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Abstract 

^^ The formation and propagation of singularities for Boltzmann equation in bounded domains has been an important 

question in numerical studies as well as in theoretical studies. Consider the nonlinear Boltzmann solution near 
Maxwellians under in-flow, diffuse, or bounce-back boundary conditions. We demonstrate that discontinuity is created 
at the non-convex part of the grazing boundary, then propagates only along the forward characteristics inside the 
domain before it hits on the boundary again. 
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; I 1 Introduction 

P^ . . . 

.^ A density of a dilute gas is governed by the Boltzmann equation 

^ dtF + vV^F^Q{F,F), (1) 

a where F{t, x, v) is a distribution function for the gas particles at a time t > 0, a position x G O C M"^ and a velocity 
w e R^. Throughout this paper, the collision operator takes the form 

^H g(Fi,F2) = B{v~u,uj)Fi{u')F2{v')dujdu- / / B{v - u,iu)Fi{u)F2{v)dujdu 

r^ = Q+(Fi,F2)-g_(Fi,F2), (2) 

^\ where u' = u + [{v — u) ■ lo]uj, v' ~ v — [{v — u) ■ uj]uj and 
^H 

• V ~~ u 

t^ B{v-u,uj) = \v-upqo{- r-w), 

O \v-u\ 

,__! with < 7 < 1 (hard potential) and 

• ^H / qo{u ■ uj)duj < +00, (angular cutoff) (3) 

^ for all u e S^ 

If the gas is contained in a bounded region or flows past a solid bodies, the Boltzmann equation must be accompanied 
by boundary conditions, which describe the interaction of the gas molecules with the solid walls. Let the domain fi be a 
smooth bounded domain. We consider three basic types of boundary conditions ( 1 1 3 ) , [2 ^ ) for F{t, x, v) at {x, v) G dflxM.^ 
with V ■ n(x) < 0, where n(x) is an outward unit normal vector at x : 

1. In-flow injection boundary condition : incoming particles are prescribed ; 

F{t,x,v)^G{t,x,v). (4) 

2. DiflFuse reflection boundary condition : incoming particles are a probability average of the outgoing particles ; 

F{t,x,v) = Cf^li{v) / F{t,x,v'){n{x) ■ v'}dv' , (5) 

J v' ■n{x)>0 

with a normalized Maxwellian ^ — e ^ , & normalized constant c^ > such that 

^i{v')\n{x)■v'\dv' = 1. (6) 

i;'-n(2:)>0 
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3. Bounce-back reflection boundary condition : incoming particles bounce back at tlie reverse tlie velocity ; 

F{t,x,v)=F{t,x,-v). (7) 

The purpose of this paper is to investigate possible formation and propagation of discontinuity for the nonlinear 
Boltzmann equation under those boundary conditions. In order to state our results, we need following definitions. 

1.1 Domain 

Throughout this paper, we assume the domain il C M'^ is open and bounded and connected. For simplicity, We assume 
that the boundary dfl is smooth, i.e. for each point xq G dft, there exists r > and a smooth function ^^o : K^ ^ K 
such that - upon relabeling and reorienting the coordinates axes if necessary - we have 

nnB{xo,r) = {xe B{xo,r) : X3 > $a;o(a;i, 2:2)}- 
The outward normal vector at dfl is given by 

1 



(8) 



n{xi,X2) 



v/l + |V,$(a:i,a:2)P 



( 9j;i$j;„(xi,a;2), 9^:2 *:eo (2:1 > 2:^2), -1 )■ 



Given {t,x,v), let [X(s),y(s)] ~ [X{s;t,x,v),V{s;t,x,v)] — [x — {t — s)v,v] be a trajectory (or a characteristics) for 
the Boltzmann equation ([l| : 

as as 

with the initial condition : [X{t; t, x, v), V{t] i, x, v)] — [x, v]. 

Definition 1 For {x,v) e J7 x M.^ , we define the backward exit time, t\y{x,v) > to be the last moment at which the 
back-time straight line [X(s;0,x,v),V{s;0,x,v)] remains in the interior of f2 : 

tb{x, v) = sup({0} U {t > : X - sv e n. for all < s < r}). 

We also define the backward exit position in dfl 

x\,{x, v) — X — tiy{x, v)v G dfl, 
and we always have v ■ n{xhix,v)) < 0. 



1.2 Discontinuity Set and Discontinuity Jump 

We denote the phase boundary in the phase space J7 x M'^ as 7 = dfl x M^, and split it into outgoing boundary 7+, the 
incoming boundary 7_ , and the grazing boundary 70 : 



7+ = {{x,v)edn 

7- = {{x,v)edQ 
7o = {{x,v)edn 

We need to study the grazing boundary 70 more carefully. 



n{x) ■ V > 0}, 
n{x) ■ V < 0}, 
n(x) ■ V = 0}. 




Figure 1: Grazing Boundary 70 
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Definition 2 We define the concave (singular) grazing boundary which is a subset of the grazing boundary 70 .' 

7f = {ix,v) e 70 : tb{x,v) ^ andti,{x,-v) ^ 0}, 

and the outward inflection grazing boundary in the grazing boundary 70 .' 

7q = {{x, w) e 7o : t\y{x, v) ^ and ti,{x, —v) — and there is S > such that x + tv E ff^ for t E (0, S)}, 

and the inward inflection grazing boundary in the grazing boundary 70 ; 

7q^ = {{x, w) G 7o : tYj{x, v) = and tb(a;, ~v) ^ and there is 6 > such that a; — rw G fj'^ for r G (0, S)}, 

and the convex grazing boundary in the grazing boundary 70 .' 

lo = {(^T^) £ 7o '■ tbix,v) — and ti,{x,—v) ~ 0}. 

It turns out that the concave (singular) grazing boundary 7q is the only part at which discontinuity can be created and 
propagates into the interior of the phase space il x M"^ . 

Definition 3 Define the discontinuity set in [0,oo) x f2 x M'^ as 

£> = |(0,oo) X [7f U7o^U7o^+ ]}u{ (i,x,t;) G (0,00) X {nxM3U7+}: t>tb{x,v)and {xb{x , v) , v) € -f^ ^ (9) 

and the continuity set in [0, 00) x fj x M"^ as 

€ ^ |{0} X OxMnu|(0,oo) X [7„U7^" ]| 

U ^ {t,x,v) e {0,00) X {nxR^Uj+} : t <tb{x,v) or (a;b(x, w), w) G 7- U 7o" |. (10) 

For bounce-back reflection boundary condition case f^ we need slightly different definitions : the bounce-back dis- 
continuity set and the bounce-back continuity set are 

Dbb = £> U Ut,x,v) £ {0,00) xQxR^ ■.t>2tb{x,v) + tb{x,-v) and {xh{x,-v),-v) EJqI, 
€bb - {{0}xnxM3|u|(o,oo)x [7_U7o^- ]} 

U < {t, X, v) G [0, 00) X {fJ X R^ U 7+} : t < ti,(x, v) or [ (xb(x, v), u) G 7- U 7q^ and t < 2t\,{x, v) + tb(a;, —v) 

or [ (xb(a;,-w),-u) G 7- U7o" and (xb(a;, w), d) G 7- U 7o" ] |, 

respectively. 

The discontinuity set D consists of two sets : The first set of ^ is the grazing boundary part 70 of J). This set 
mainly consists of the phase boundary where the backward exit time ti,{x, v) is not continuous (Lemma l2|. The second 
set of ([9| is mainly the interior phase space part of 5D, i.e. 2) n {[0,oo) x fl x R'^}, which is a subset of a union of all 
forward trajectories in the phase space emanating from 7q . Notice that S does not include the forward trajectories 
emanating from ^q U 7g^ because those forward trajectories are not in the interior phase space [0, 00) x fJ x M'^. We also 
exclude the case t < ti,{x, v) from S). In fact, considering the pure transport equation, t < ib(2;, v) implies the transport 
solution at {t,x,v) equals the initial data at (a; — tv,v) and if the initial data is continuous, we expect the transport 
solution is continuous around (t, x, v). Notice that we exclude the initial plane {0} x fj x M'^ from S because we assume 
that the Boltzmann solution is continuous at i = . 

The continuity set (^ consists of points either emanating from the initial plane or from 7_ U 7o , but not 7q . 

Furthermore we define a set including the grazing boundary 70 and all forward trajectories emanating from the whole 
grazing boundary 79. 

Definition 4 The grazing set is defined as 

& = {{x,v) eClxR^ ■.n{x):,{x,v)) -v = n{x-ti,{x,v)v) -v ^0}, (11) 

and the grazing section is 

&a,^{v£R^ : {x, v) e&} = {v G R^ ; n{xi,{x, v)) ■ v ^ 0}. 
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Obviously the grazing set © includes the discontinuity set J). In order to study the continuity property of the Boltzmann 
solution we define : 

Definition 5 For a function (l){t,x,v) defined on [0, oo) x {J7 x ]R'^\©} we define the discontinuity jump in space 
and velocity 

[(t){t%,v = lim sup \(t>it,x\v')-(f>{t,x",v")\, 

^^^ (x' ,v'),{x'\v")e{nxW=^\&}n{B({x,v)i5}\{x,v)} 

and the discontinuity jump in time and space and velocity 



lim sup 

^^^ t',t" e B{t;S) 

{x',v'), {x",v") e {fix E^\0} n {-B((x, v); S)\{x, v)} 



m',x',v')-<j>{t",x",v'% 



where @ is defined in Definition ul We say a function cf) is discontinuous in space and velocity (in time and space and 
velocity) at (t, x, v) if and only if [(f>(t)]x^v 7^ ([0]t,a;,i; 7^ 0) and continuous in space and velocity (in time and space and 
velocity) at {t,x,v) if and only if [4>{t)]x,v = {[(l>]t,x.,v = 0). 

Notice that the function (jj is only defined away from the grazing set ©. Once the discontinuity jump of given function 
ip is zero at (i, x, v) then the function ip can be extended to [0, 00) x J7 x K'^ near (t, x, v). Because of those definitions 
we can consider a function which has a removable discontinuity as a continuous function. And a non-zero discontinuity 
jump [</>]t,2:,i, j^ means (j) has a "real" discontinuity which is not removable. 



1.3 Main Result 

Now we are ready to state the main theorems of this paper. In order to state theorems in the unified way we use a 
weight function 

wiv) = {l + p^\vff, (12) 

such that w^'^{v){l + \v\y e L^. 

Theorem 1 (Formation of Discontinuity) Let fl be an open subset ofM.^ with a smooth boundary dft. Assume H. 
is non-convex, 7g 7^ 0. Choose any (xo,wo) G 7o ''^*^^ ''^o 7^ 0. For any small S > 0, 

1. There exist ig G (0,min{(5, ib(a;o, — I'd)}) '^'^d an initial datum Fo(x,v) which is continuous on ft x M.^ U {7_ U 7^ }, 
and an in-flow boundary datum G{t,x,v) which is continuous on [0, 00) x {7^ U7f}, satisfying 



Fo{x,v) = G{0,x,v) for (x, w) G 7-^ U 7^, 



and 



Fo-P 



//i 



L°°(nxK3) 



sup 

tG[0,oo) 



G{t)-p 



v^ 



i°°(7-) 



<<5, 



(13) 



(14) 



such that if F on [0,cxd) x f2 x M"^ is Boltzmann solution of Uy with the in-flow boundary condition (Mj then F is 

discontinuous in space and velocity at (to,xo,vo), i.e. [F{to)]xo,vo 7^ 0. 

2. There exist to € {Q,m.in{6,tb{xo, —vq)}) and an initial datum Fo{x,v) which is continuous on ft x R^ U {j- U 79}, 

satisfying 



Fa{x,v) = c^p{v) Fo{x,v'){n{x) -v'jdv' for (a:, w) G 7_ U 79 . 

^t;'-n(a;)>0 



(15) 



w 



A* 



a/M 



L°°(nxR3) 



<S, 



(16) 



such that if F on [0, 00) xilxM.^ is Boltzmann solution of uy with the diffuse boundary condition ([5|) with the compatibility 
condition (15) then F is discontinuous in space and velocity at {to,xo,vo), i.e. [F{to)]x„,vo 7^ 0- 



3. There exist to G (0 , m in{(5, tb(xo, —vq), t^lxp, vq)}) and an initial datum Fq{x,v) which is continuous on D, x M' U 
{7_ Uj^}, satisfying [1" 



Fo{x,v) ^ Faix,~v) for (a;, w) G 7- U 73 . 



(17) 



such that if F on [0,cx)) x J7 x M^ is Boltzmann solution of tw with the bounce-back boundary condition rTp then F is 
discontinuous in space and velocity at (to,Xo,Vo), i.e. [F{to)]xo:va ¥" 0- 



The smallness of given data (14), (16) ensures the global existence of Boltmzann solution for all boundary conditions 



[13j . Notice that we can observe the formation of discontinuity for any point of the concave (singular) grazing boundary 
7^ of any generic non-convex domain il. If we assume that Fq{x, v) is conti nuo us on il x ffi^U {7- U 7q } and G{t, x, v) 
is continuous on [0, 00) x {7. U 7q }, and that the compatibility conditions (13), (15) and (17) are valid up to 7_ U 7g , 
then Theorem 1 implies the continuity breaks down at the concave(singular) grazing boundary 7q after a short time 
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to € {0,m.m{6,tb{xo,—vo)}) for in-flow Q, diffuse ([s]) boundary condition and to € {0,mm{S,tb{xo,—vo),tb{xo,vo)}) 
for bounce-back ^ boundary condition. For this generic cases, we said the Bohzmann solution F has a local-in-time 
formation of discontinuity at {to,Xo,vo)- 

Once we have the formation of discontinuity at {to,XQ,vo) G 7o, we further estabUsh that the discontinuity propa- 
gates along the forward characteristics. 

Theorem 2 (Propagation of Discontinuity) Let fl be an open bounded subset ofM.^ with a smooth boundary dfl. 
Let F{t,x,v) be the Boltzmann solution of Iw with the initial datum, Fo which is continuous ok il x M"^ U {7- U 70}? 
and with one of the following boundary co nditi ons : 

1. For in-flow boundary condition M, let Ujw and {I4) be valid and G{t,x,v) be continuous on [0, cx)) x {7_ U7q}. 

2. For diffuse boundary condition (m, assume (16) and (J5[ ). 

3. For bounce-back boundary condition O), assume \16\ and [F, 
Then for all t G [to, to -\- ti,{xo, —vo)) we have 



where Ci > only depends on 



fFh , /. . ^ < p-Ci(i+|fo|)"'(t-to)rp,'+ )i 



(18) 



/7I ML°°([0,oo)xnxM3)' 

On the other hand, assume [F{to)\xa,vo 7^ 0, and to £ (0, tb(xo, — Wq)) for in-flow and diffuse boundary conditions 
and to e (0,Tam{t\3{xo, —vo),t\:,{xo,vo)}) for bounce-back boundary condition, and a strict concavity of dfl at xo along 
vo, i.e. ^^ 

(19) 



i,3 



{vo)idx,d^i'^ixo){vo)j < -C^ 



Then for all i G [to , ^0 + 
{t,xo + {t-to)vo,vo), i.e. 



th{xo,^vo)), the Boltzmann solution F is discontinuous in time and space and velocity at 
[F]t,xo+(t-to)va,vo 7^ and 



where < C < 1, and C2 — C^U 



w ~^ r 00 ) £ ^ which is positive for sufficiently small 



(20) 



yTJi Ml°°([0,oo)xS1xK3)' 



The strict concavity condition ( 19 1 rules out some technical issue of the backward exit time tb- Our theorem characterize 



the propagation of discontinuity before the forward trajectory reaches the boundary. In the case that the forward 
trajectory reaches the boundary, i.e. t > io + ^bl^^oi ""^^o); the situation is much more complicated. Denote xi = 
2^0 + ih{xoT—vo)vo,ti — to -\- t-t,{xo,—vo). If the trajectory hits on the boundary non-tangentially, i.e. (a;i,Uo) € 7+, 
for in-flow and diffuse boundary cases, the discontinuity disappears because of the continuity of the in-flow datum and 
the average property of diffuse boundary operator. For bounce-back case the discontinuity is reflected and continues 
to propagate along the trajectory. If the trajectory hits on the boundary tangentially, i.e. (a;i,wo) G 70, there are 
three possibilities. Firstly if (a;i,wo) G 79^ then the situation is same as the case (a;i,vo) G 7+ above. Secondly, if the 
trajectory is contained in the boundary for awhile, i.e. there exists 5 > so that Xi -\- svo S dfl for s E (0, S) then it is 
hard to predict the propagation of discontinuity in general. Assuming certain condition on fl for example Definition |6l 
we can rule such a unlikely case. 

The last case is that (a;i,wo) G 7o- Assume we have a sequence of {t„ — tn-i + ti,{xn-i,—vo)} and {x„ = 
Xn-i -\-ti,{xn-i, — wo)'^o} G 90 so that {xn,vo) G 7^ , and a directional strict concavity (19 1 is valid for each (x„,tio). We 
can show the propagation of discontinuity also between the first and the second intersections, i.e. [-F]t,a;o(t-to)i'o-i'o 7^ '-' 
for t G [^1,^2) in general. For t > t2, if we have very simple geometry, for example the first picture of Figure 2, we can 
show the propagation of discontinuity, i.e. [F]t^xo(t-to)vo,vo 7^ for i G [tn,tn+i) even for n — 2,3. But in general, for 
example the second picture of Figure 2, we cannot show [F]f,ca{t-ta)vo.va 7^ for t G [tn,tn+i) for n>2. 





Figure 2: Grazing Again 

The next result states that Theorem [I] and Theorem [2] capture all possible singularities (discontinuities), despite non- 
linearity in the Boltzmann equation. In other words, the singularity of the Boltzmann solution is propagating as the 
linear Boltzmann equation and no new singularities created from the nonlinearity of Boltzmann equation. 
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Theorem 3 (Continuity away from S) Let fl be an open bounded subset o/M"^ with a smooth boundary 9J7. Let 
F{t,x,v) be a Boltzmann solution of (Mj with the initial datum Fq which is continuous on fl x M.^ U {7- U 7+ U 7o~} 
and with one of 
1. In-flow boundary condition (Op. Assume (I4) is valid and the compatibility condition 



Fo{x,v) = G{Q,x,v) for (a;, f) G 7_ U 7^-, (21) 



and G{t,x,v) is continuous on [0, cx)) x {j_ U 7q }. 

2. Diffuse boundary condition [^. Assume Um is valid and the compatibility condition 



Fo{x,v) ^ Cfj,fi{v) Fo{x,v){n{x) -v'jdv' for {x,v) e'-f-Dj^ . (22) 

J v' ■n{x)>0 



3. Bounce-back boundary condition (MJ. Assume (16) is valid and the compatibility condition 



Fo{x,v)^Foix,-v) for {x,v) e j- U^o' . (23) 

Then F{t,x,v) is a continuous function on C for 1,2 and a continuous function on €bb for 3. If the domain Vt does not 
include a line segment (DefinitionW^ then the continuity set £ and £;,(, are the complementary ofT) and T)bb respectively. 
Therefore F(t,x,v) is continuous on (S)"^ for 1,2 and continuous on {TibbY f'^^ ^■ 

Definition 6 Assume f2 e M"^ be open and the boundary dVt be smooth. We say the boundary dH- does not include a 
line segment if and only if for each Xq G dfl and for all (ui, 1*2) G S^ there is no S > such that 

^Xo{tUi,TU2) 

is a linear function for t G {~S, 6) where ^xq from H). 

1.4 Previous Works and Significance of This Work 

There are many references for the mathematical study of different aspects of the boundary value problem of the Boltz- 
mann equation, for example [in][IS][IS] and the references therein. In [T^], an unified L^ — L°° theory in the near 
Maxwellian regime is developed to establish the existence, uniqueness and exponential decay toward a Maxwellian, for 
all four basic types of the boundary conditions and rather general domains. 

The qualitative study of the particle-boundary interaction in a bounded domain and its effects on the global dynamics 
is a fundamental problem in the Boltzmann theory. One of challenging questions is the regularity theory of kinetic 
equations in bounded domain. This problem is hard because even for simplest kinetic equations with the differential 
operator v ■ S/^, the phase boundary dfl x M.^ is always characteristic but not uniformly characteristic at the grazing 
set ^Q — {{x,v) : X £ dfl, and v ■ n(x) — 0}. In a convex domain a continuity of the Boltzmann solution away from 
7o is established in [13J for all four basic boundary conditions. In a convex domains, backward trajectories starting at 
interior points of the phase space cannot reach points of the grazing boundary 70, due to Velocity Lemmaf |llj | 15 | ). 
where possible singularities may exist. 

In general, on the other hand, in a non-convex domain, backward trajectories starting at interior points of the phase 
space can reach the grazing boundary. Therefore we expect singularities will be created at some part of grazing boundary 
7o and propagate inside of the phase space. This question has been attracting a lot of attentions from early '90s, see 
references in pp. 91-92 in Sone's book ^50". For Boltzmann equation, most of works are numerical studies pD] | 21 ) P^ 
and few mathematical studies. 

Once we enlarge our survey to propagation of singularities which already exist on initial data or boundary data, 
there are some mathematical works [3] [5] [6] [7] [8] as well as numerical works [7j20]. In [3 , for linear BGK model, a 
propagation of discontinuity ,which exists already in the boundary data, is studied mathematically and also numerically. 
In [S], for the full Boltzmann equation in the near vacuum regime, a propagation of Sobolev H^/"^^ singularity, which 
exists already in the initial data, is studied and same effect has been recently shown in the near Maxwellian regime 

In Vlasov theory, we refer to (2] [9] [25] for the boundary value problem. Singular solutions were studied in [11 exten- 
sively. In '11] , the non-convexity condition of boundary is replaced by the inward electric field which has a similar effect 
with non-convexity of the boundary. In convex domains, Holder estimates of Vlasov solution with specular reflection 
boundary is solved recently |15| | 16 ) . but Sovolev-type estimate is still widely open. 

Our results give a rather complete characterization of formation and propagation of singularity for the nonlinear Boltz- 
mann equation near Maxwellian in general domain under in-flow, diffuse, bounce-back boundary conditions. There is 
no restriction of the time interval. More precisely we show that for any non-convex point x of the boundary and velocity 
tangent to 90 at x, there exists an initial datum (and in-flow datum, for in-flow boundary condition case) such that 
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the Boltzmann solution has a jump discontinuity at {x,v). Once the discontinuity occurs at the grazing boundary, this 
discontinuity propagates inside along the forward trajectory until it hits the boundary again. And except those points, 
the grazing boundary and forward trajectories emanating from the grazing boundary, we can show that the Boltzmann 
solution is continuous. (Continuity away from 23) 

1.5 Main Ingredients of the Proofs 

1. The Equality induced by Non-Convex Domain 



We consider near Maxwellian regime and linearized Boltzmann equation ( 32 1 . The formation of discontinuity is a 
consequence of following estimate. Assume (x, v) G 7g as below picture so that for sufhciently small i > the backward 
trajectory x — tv is in an interior of the phase space. For simplicity we impose the trivial in-flow boundary condition 



G{t,x,v) = ^jiiv) which corresponds g{t,x,v) = (92). Consider points (xj^,f^) in 7_ and {x'^,v'^) missing the non- 



convex part near (x, v) and both sequences converge (x, u) as ti — > oo. 




Figure 3: 

Now suppose the solution / of the linearized Boltzmann equation be continuous around {x,v). Then the Boltzmann 
solution / at {x'.^,v'^) 

and at (a;^,u"), 

Jo 
converges each other as n -^ cxd. Then we have the following equality 

fo{x -tv,v) = - f e''(")^{X/ + r(/, /)}(5, x~{t~ s)v, v)ds. (24) 

Jq 

Thanks to |13j, the pointwise estimate of /, with some standard estimates of K, T, the right hand side of above equality 



has magnitude O(t)||/o||oo(l + ||/o||oo)- If you choose fo{x — tv,v) — ||/o||oo then the above equality (24) cannot be true 
for sufficiently small t unless the trivial case /o = 0(F = /x). Therefore the Boltzmann solution / cannot be continuous at 
{x, v). For diffuse ([5]), bounce-back ^ boundary conditions we also obtain the equality induced by non-convex domain 



similar as (24) 



This argument bases on the idea that free transport effect is dominant to collision effect if time f > and the per- 
turbation -^^ is small. 

2. Continuity of the Gain Term Q+ 

The smoothing effect of the gain term Q^ is one of the fundamental features of the Boltzmann theory. There are lots 
of results about the smoothing effect in Sobolev regularity, for example 

iig+(0,v)ii^«^ <cm\LAmL^, 



with some assumption on various collision kernels |18j |26j |27j . To study the propagation of singularity and regularity, in 
the case of angular cutoff kernel (l3]), it is standard to use Duhamel formulas and combine the Velocity average lemma 
and the regularity of Q+ ^\. For detail, see Villani's note [23] especially pp. 77-79. 

In order to study the propagation of discontinuity and continuity we need a totally different smooth effect of Q+. 
For the discontinuity induced by the non-convex domain, we need following : Recall the grazing set 25 in Definition |4] 
A test function (j)(t,x,v) is continuous on [0,r] x (fi x K^)\e5 and bounded on [0,T] x ft x K^. Then 

Q+i<l),(l)){t,x,v) e C°([0,r] xrJxM^). (25) 
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Recall that the grazing set © = {{x, v) eQxR:v e (Sx}. The grazing section &j; ^ {ru G M^ : i > 0, m G ©a; n S^} is 
a union of straight lines in velocity space M"^ and two dimensional Lebesque measure of <&x H S^ is zero (Hongjie Dong's 
Lemma, Lemma 17 of |13j). Moreover, using continuous behavior of C^. in a;, one can invent a very effective covering 
of &x (Guo's covering. Lemma 18 of [13]). Because of those geometric and size restriction on ©, even the gain term 
(3+ is an integration operator in v alone, we can prove the smoothing effect of Q+ in C°([0,T] x i7 x W") for t^x and 
V, see Theorem Wl Notice that those smoothing effect on C^^,,, has been beheved to be true for long time without a 
mathematical proof in numerical communities [1], pl587 of [3], p502 of [2T] . 

The main idea to prove the smoothing effect in C^^^ is to use the Carleman's representation for Q^((j),(j))(t,x,v) 
which has been a very effective tool [H] [55] [37] . 

0(^,a;,«') ^_^ / 4'{t,x^v[)B{2v-v' -v[,——^)dv[dv', (26) 

with the hyperplane Ey^i — {v'l G M.^ : {v[ — v) ■ {v' — v) =0}. We will show the smallness of 

\Q+{<P, <P){t, X, v) - Q+{<j>, (p){t, x,v)\, 

for \{t,x,v) — {i,x,v)\ < 6. Assume we have sufficient decay of 4> for large v. Replace the integrable kernel , _^^,|2 

by smooth compactly supported function and cut off the singular part of B(2v — v' ~ v[, Zi_yn ) to control the above 
difference as 

0{S)M\lo + C f m,x,v')-cf>{i,x,v")\ f \<P{i,x,v'l)\dv'ldv' 



\v'\<N JE^^,,r]{\v'^\<N} 



+ C \(j){t,x,v')\{ <j){t,x,v[)dv[^ (j){t,x,v';)dv'l\dv', 

J\v'\<N [J E^^,n{\v{\<N} •^-EB„"n{|i';'|<^'} J 

where v"{v') is chosen to be v' — (v ~ v) for convenience. 

One can easily control the integration at the first line. Because for the first term, integrating over w', we can cut off 
a small neighborhood of (Sx from \v'\ < N. Away from that neighborhood, using the continuity of away from 0^ we 
can control the integrand pointwisely. 

In order to control the second line integration we have to control the difference in big braces. To do that we choose 



a special change of variables for w" ,(40). Under this change of variables the second line is bounded by 



C \(f){t,x,v')\ \<P{t,x,v[)-<p(t,x,v'l)\dv[dv'. 

J\v'\<N JE^^.,n{\v[\<N} 

The second integration term above is a function of t,x,i,x and v. Unfortunately one cannot expect a pointwise control 
(smallness) of the second integration for all v : even the grazing section ©3., where (f>(t,x,v'i) might have discontinuity, 
is small, i.e. 2-dimensional Lebesque measure of ©a, n S^ is zero, the measure on the plane Eyy' could be large(even 
infinite). However, in Section 3.3, we can show that that bad situation happens for very rare v' in {v' E M.^ : \v'\ < N} 
and use the integration over v' to control the above integration. 

3. New Proof of Continuity of Boltzmann solution with Diffuse Boundary Condition 

In Section 5.2 we prove a continuity away from S of Boltzmann solution with diffuse boundary condition using simple 



iteration scheme (1031 with iteration diffuse boundary condition (131). This iteration scheme has several advantages. 
First it preserves a continuity away from 2) as m increasing, that is, if /i™ is continuous away from 1) then /i'"+^ is also 
continuous away from S). Second, the sequence {/i™} has uniform L°° bound and moreover it is Cauchy in L°° for in-flow 
boundary condition /i'"!^. = wg. Therefore h — lim/i'", a solution of the linear Boltzmann equation is continuous local 
in time. Combining with uniform-in-time boundedness of Boltzmann solution f |13|). we achieve the continuity for all 
time. In order to apply this idea to diffuse boundary condition, we use Guo's idea [13j : A norm of the diffuse boundary 
operator is less than 1 effectively, if we trace back several bounces. This approach gives simpler proof for the continuity 
of Boltzmann equation with diffuse boundary condition with convex domain (see Lemma 23 25 of |13j). 

1.6 Structure of Paper 

In Section 2, we state some preliminary facts which are useful tools for this paper. In Section 3, we state and prove 
the continuity of Q+ (Theorem 4). In Section 4 6, we deal with in-flow boundary, diffuse boundary and bounce-back 
boundary, respectively. For each section, first we prove the formation of discontinuity (Theorem 1). Then we show the 
continuity away from S) (Theorem 3). Using this continuity, combining with continuity of Q^, we show the propagation 
of discontinuity (Theorem 2). 
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2 Preliminary 

In this section we study continuity properties of the backward exit time tb{x,v) and, a measure theoretic property and 
geometric covering of the grazing set ©, and estimates of Bohzniann operators and the Carleman's representation. 

We use Lemma 1 of [13 , basic properties of the backward exist time tb{x,v) : 

Lemma 1 \13f Let fl be an open bounded subset of M.^ with a smooth boundary dfl. Let {t,x,v) be connected with 
{t — tb{x, v), Xb{x, v), v) backward in time through a trajectory of (UTTI). 
L The backward exit time ti,{x,v) is lower semicontinuous. 
2 If 

vn{xb{x,v)) < 0, (27) 

then {t\3{x,v),Xb{x,v)) are smooth functions of {x,v) so that 

Va^tb — 7 T, Vt,tb — -. r, 

V ■ n(Xi,) V ■ n(x-h) 

VajXb = / + V^^ib <8) V, Vi,Xb = thi + V„ib ® V- 



For a convex domain, if a point (x, v) is in the interior of the phase space, i.e. (a;, w) G 51 x M'^ then the condition ( 27 ) is 
always satisfied and hence tb(a;, v) is smooth due to Lemma fl] However for a non-convex domain, there is a point {x, v) 
in f7 X M'^ but (a;b(a;, v,),v) G Jo, i-e. v ■ n{x]3{x, v)) = 0. We further investigate a continuity property of tb for that case. 
Indeed, discontinuity behavior oit\,{x,v) for {xh{x,v),v) G 7g is a main ingredient of the formation of discontinuity. 

Lemma 2 Let fl ^ M.^ be an open set with a smooth boundary dfl. Assume (a;o,uo) G fJ x M'' with uq 7^ and 
tbixo,vo) < oo. Consider (xqj'Uo) e 6, i.e. (a;b(xo, wq), wo) G 7o- 

If ixh{xo,Vo),VQ) £ Jq^ thenth{x,v) is continuous around {xq^Vq). 
If (a;b(a;o, wo); ^o) G 7o then ib(a^, v) is not continuous around (xq, vq). 

Recall J Q^ and 'Jq in Definition^ 

Proof. Throughout this proof, without loss of generahty we assume that dO. is a graph of $ locally and $(0, 0) = and 
(i92;i$,9x2^)(0, 0) — (0,0). Moreover assume Xq = (|a;o|, 0,0),t'o = (|uo|:0,0) and ib(a;oji'o) = ]f4 ^^ ^^^^ x\3{xo,vo) = 
(0,0,0) = (0,0, $(0,0)). 

First, let ixb{xo,VQ),vo) G 7q .By the definition of 7q , we have $(— t, 0) > and ^(t, 0) < for < t << 1. 
Using the continuity of <i>, choose sufficiently small e > 0, 6 > such that $(—(5, y) > | and ^{6, y) < — | for < |y| < S. 
Fix X = (xi, 2:2,2:3) '^ xq and v — {vi,V2,V3) r^ vq. We define 

^(x, V, t) — x^ — tv^ — ^{xi — tVi, X2 — tV2). 

FOI t' = ^, ^'ix,V,t'o) ^ ~^S,X2 - ^V2) + X3 - ^V3 > '^ ior ixi,X2,X3) ^ {\xo\,0,0), {VI,V2,V3) ^ {\vo\,0,0). 

For t" = ^, '${x,v,t") = -^-S,X2 - ^V2) + X-, - ^V3 < -f for {x,,X2,X3) ~ (|a;o|,0,0), {v,,V2,V3) ~ 
(|wo|,0,0). Using the continuity of $ and ^P, there exists ^* G (f^- — ^, f^- + :^) so that ^{x,v,t^) = 0, i.e. tb{x,v) — i*. 
If j: - xo and u - uo then 2^ - ^ - |^ - r^ = ib(a;o, wo) - A and 2^ + ^ - 1^ + A ^ tbixo, vq) + A so that 
i» G {tbixo,Vo)- 4;|,ib(a:o,Wo) + p^y)- 

Next, let {xb{x,v),v) G 79 ■ By the definition of the concave grazing boundary Jq, we have $(— r, 0) > and 
i'(T,0) < for < T << 1. Choose a sequence x„ — (|a:o|,0, ^). There exists e > such that tb(a^n,wo) > tb{xo,vo) +£ 
for sufficiently large n. This implies that (a:„,wo) —?> (a;o,wo) but ib(2:n,wo) ^^ ti,{xQ,vo) as n — > 00. ■ 

In the next two lemmas, we consider the grazing set © (Definitional including the discontinuity set S). The next 
lemma. Lemma 17 of [T3j due to Hongjie Dong, is important to control a size of &. We denote m2 as a standard 
2-dimensional Lebesque measure and TO3 as a standard 3-dimensional Lebesque measure. Recall that the grazing section 
&x in Definition |4] 

Lemma 3 Ucif If dCl is C^ then the grazing section ©j. restricted to §^ has zero 2-dimensional Lebesque measure, i.e. 

m2(©,n§2)=0, 

for all X ^ f2. 

With condition 7712 (©^^ fl §^) =0, we can construct the Guo's covering which is little bit stronger than the original one 
in Lemma 18 in [13] . 
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Lemma 4 (Guo's covering) fl^ Assume m2{<&x n §^) = is valid for all a; e fj. Let Bn = {v eM? : \v\ < N}. 
Then for any e > and N.^, > there exist Se^N,N, > 0, and le,N,N,,n balls B{xi; ri), B{x2', r2)---, B{xi; ri) C Vt, as well 
as open sets Oxi,Ox2,---Oxi of B^ which are radial symmetric, i.e. 

O^^ = {iu e M'"^ : t > 0, ve O:,. n S^}, 

with m^{Oxi) < -^ and m2{Oxi H S^) < j^2j^ for all 1 < i < le.N,N,,ni such that for any x G 0,, there exists Xi so that 
X € B{xi; ri) and for v ^ Ox^ , 

\v ■ n[xb{x,v))\ > 6^^N,N, > 0, 

or equivalently 

Ox,^ IJ {v(^Bn ■.\v-n{xb{x,v))\<6e,N,N.}^ IJ ©a^nSjv. 

x^B{xi\ri) x^B{xi:ri) 

Combining Lemma [3] and Lemma |4j we have following lemma, which is useful to prove Theorem [4j Namely, a function 
which is continuous away from the grazing set is uniformly continuous except arbitrary small open set containing 0. 

Lemma 5 Assume (j) is continuous on [0,T] x (il x {w G M'^ : -p < |w| < Af})\©. For fixed x E fl and e > and 
N^ > 0, there exist 

d = 6{<j>,n,s,N,,x,j^,N)>0, (28) 

and an open set Ux C {v Q M.^ : jj < \v\ < N} which is radial symmetric, i.e. Ux — {tv eM"^:f>0, v E Uxf^ S'^} 
with m^{Ux) < -j^ and m2{Ux n S'^) < j^^^ such that 

- s 

\(t){t,x,v) -(/)(i,x,-u)| < — , 

forv e {weM^ : ^ < \v\ < N}\Ux and \{t,x,v)- {t,x,v)\ < S. 

Proof. Let x ^ x. Due to Guo's covering [T^], LemmaWl we can choose B{xi; ri) including a; and x, as well as Oxi C M^ 
so that 

Ox^^ U ©y n Bjv D IJ (^yH Bn, 

yeB(xi-,ri) yeB(x;5) 

with 7713 (OxJ < -^. Notice that mz[Oxi) — m^iOxi). We can choose an open set Ux^ so that mz[Uxi) < 2m3(Oa;J and 
Oxi C Uxi- Since both of Ox^ and B^XUx^ are compact subsets of -Bat, we have a positive distance between two sets, 
i.e. 

= inf{|C - ei : C e O,, and ^ e Bn\UxA- 

Assume 6 < 3/2. Fix x e (} and v e {v e M.^ : jj < \v\ < N}\Ux. Then \{x,v) - {x,v)\ < 6 implies that 
V & {v G M.^ : jj < \v\ < N}\Oxi. For such x,v,x and v consider the function (j) as it's restriction on a compact set 
[0,T] X B{x;S) X Bi^\Oxi. Therefore 4>[id,T]y.B(x:S)y.BN\o^. i^ uniformly continuous. Hence \(j){t,x,v) — 0(t,a;,w)| can be 
controlled small uniformly if (5 > is sufficiently small. ■ 

We will use the Carleman's representation [12] [26] in the proof of Theorem [4] crucially. Let Q+{(j),ip) be defined by 
([2]) and let ip = ^^{v) and (f) = (t>{v), v €W^ make Q+{'ip, (jj) < 00 almost everywhere. Then the Carleman's represen- 
tation is 

g+(V, mv) = 2 / ^{v') \ / cf>{v[)B{2v -v'-v[, ^^^y^)dv[dv', 



(29) 



where E^^i is a hyperplane containing w G M'^ and perpendicular to Z,J"^i G §^, i.e. 

£:^„, = {v[ e M^ : {v[ - v) ■ {v' -v)= 0}. (30) 



In the proof of Theorem [4] we need to control the integration over Eyyi in ( 29 ) frequently : 
Lemma 6 For a rapidly decreasing function : M-|_ — > M-|_, we have 

<l>{W\)B{2v -v'- V,', ^^, ~ l^[,^ )dvi' < Q(l + \v- v'D, (31) 

where C^f, only depends on (p. 



3 CONTINUITY OF THE COLLISION OPERATORS 11 

Proof. For fixed v' and v, let us denote {61,62,63}, with §3 = Z,Zy\ , be the orthonormal basis of M."^ such that any 



v[ & Ey^i can be written as Uj = u + 77161 +7/262- Since u' — u _L £'„„' from (30), there is 773 such that 77' — 77 = 7/363 where 
l^l = It' — w'l- Then we can write 2v — v' — v[ — v — v' + v — v'i = —77161 — 77262—77363 and \lv — v' — v'-^ = 7yi+77| + |v' — t;^. 



Moreover 7;' — Wi = —ri\e.\ — 7/262 + 7/363. We can write the left hand side of (31 1 as 



00 /•CO 



<^(^?+^?2 + bn 



< 



00 J —00 

00 /"OO 



-V3 



-m 



7/^ + 77I + |w — i;'p 
<^(^? + vlM + r,l + \v' - v\^)i-\r,l + 77I - \v' - v\')d7j^dm 



-7/2 I ^771^772 



-00 >/ —00 

00 /'OO 



J ~oo J —00 

< C^il + lv'-vP). 
■ In terms of the standard perturbation / such that F = 11 + ^/Jlf, the Boltzmann equation can be rewritten as 

{dt+v-V + L}f^ r(/, /), /(O, X, v) = /o(x, v), 
where the standard linear Boltzmann operator, see J12j . is given by 

L! = uf- Kf, 
with the collision frequency v{v) = J jw — u\'' ^{u)qo{ |"I"| • uj)dujdu for < 7 < 1 and 

1 



C, 



-(i + |i'i)^<K«)<a(i + kir, 



if/ 

r(/,/) 



k{v,v')f{v')dv' 



^ 0+(^,V/I/) + ^Q+(y^/,^)-^Q_(y^/,^), 



//i 



//x 



/H 






We recall two estimates of operators K and T from 13J. The Grad estimate for hard potentials : 



|k(w,w')l <C^{\v-v'\ + \v-v'\-^}e' 



H"-^'! 



„'|2|2 



i„-vi- 



Recall w in (12). Let < 6* < |. Then there exists < e{9) < 1 and Cg > such that for < e < e{9), 
{|t; - v'\ + \v~- v'\-^}e'~^ ' "^ i--"'i^ \^ .,,,„ d7;' < 



Ck 



7«(7;')e«l^'l' ~ l + |f| 



For the nonlinear collision operator 



Also we recall a standard estimate 



for 6e L^ 



\wTigi,g2){v)\ < Cr(l + \v\r\\wgi\\oo\\wg2\\c 



{v')\v-v'\'-'dv' -. {l + \v\y, 



3 Continuity of the Collision Operators 

In this section we mainly prove the following Theorem : 

Theor6m 4 (Continuity of Q+) Assume F{t,x,v) is continuous on [0,T] x (fl x M.'^)\<8 and 

l|w)"^i^|L~([0,T]xnxR3) < +00, 

where w = ^f 2\l\2)i3 with p G M and /3 > 0. Then Q+iF, F){t, x, v) is continuous in [0, T] x 51 x M'^ and 

sup \i^^^'w^^Q+{F,F){t,x,v)\ < 00. 

[0,T]xnxK3 



(32) 



(33) 

(34) 
(35) 



(36) 
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Theorem Wl a smooth effect in C°^ ^, is the crucial ingredient to prove Theorem and Theorem pi This smooth effect 
of the gain term ensures that there is no singularity created by the nonlinearity of Botlzmann equation. 



Proof of (36 1 . It is easy to show the boundedness (36) from 
iy~^w~^Q+iF,F){t,x,v) < 



i'{v)w{v) V Js 



B{v — u,uj)w{u')'w{v')dujdu X \\w ^F\ 



< i^ivy^ 



B{v - u,uj)—- — g 

g2 {1+P^\U\^)P 



dtddu X \\w ^F\ 



< C i'{v)-^iy{v)\\w''^F\\l^<C\\w-'^F\ 



L~([0,T]x(nxR3)) 



where we used (35) and \u' 



Next wc will show the continuity part of Theorem |4] The goal of following three subsections is to show 

For fixed e > and {t, x, v) £ [0,T] x il xR^ , there is (5 > such that 
\Q-^-{wh,wh){t,x,v) — Q^{wh,wh){t,x,v)\ < e for \(t,x,v) — {t,x,v)\ < S. 



(37) 



3.1 Decomposition and Change of Variables 

In this section, we use the Carleman's representation to split Q^{iuh, idh)(t, x, v) — Q^{il]h, wh){t, x, v) in a natural way 
(38), and introduce two change of variables (39) and (40). 



It is convenient to define 



h = w^^F, 



where 1 1 /i lino = 



L°°([0,T]x(axR3)) 



Representation ([29| we have 

Q+{wh, wh){t, X, v) — Q+(wh, wh)(t, x, v) 



^F\\L'=°([o,T]x(nx'R^))- Choose {t,x,v) ^ {t,x,v). Using the Carleman's 



iu{v")h{t,x,v") 



," 



^^ / w{vT)h(i, X, v'l)B{2v - v" - vl ^^jtAt,) dv'Idv" 
\v - v"Y Je- „ \v - "'^iV 



■V — 

B 



w{v')h{t, X, v') 



\v - u'|2 



w{v[)h(t, X, v[)B(2v — V — u'l , ^— tt) dv\dv 

\v' — v\\ 



2/ {A- A} [ Bdv'{dv" + 2[ A' [ {B-B'} dv'Idv'. 



(38) 



In order to control the first term of (38), we need to compare arguments of A and A! . For that purpose, we introduce 
the following change of variables : 



Lemma 7 For fixed v and v in M? we define 

v" = v" {v'] V, v) = v' — {v ~ v). 
Then two planes Eyy" and Eyyi have same normal direction. The distance between to planes is \{v — v) ■ i","", | 



(39) 



Proof. Assume (|39|) . Clearly "g^, ' ~ I where / is 3 x 3 identity matrix. The normal direction of i^c-u" is |^„_^| = |^,_^| 



which is also the normal direction of E^yi . To measure a distance between two planes E^yi and Eyyii, we consider 
the line passing v and directing ,^,Zy\ 1 which is v{s) — |^,~"| S + v. The solution of f (s*) G Eyy" is a solution of 

= {v" — v) ■ {v{s) — v) — {v' — v) ■ {v{s) —v) = \v' — v\s + [v' — v) • {v — v). Easily we have s* = 7i/-d~^ • ^i'^'^^ ^(^) 
is unit-speed line we know that \v{s^) — v{0)\ is the distance between Eyy'^ and Eyy/. m 



An important property of ( 39 ) is that two planes Eyy/i and Eyyi have the same normal direction. In order to con- 



trol the second term of (38), we need to compare arguments of B and B' especially v'l G Eyyi and v'{ e Eyyii. For that 



purpose, we introduce the following change of variables 
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Lemma 8 For fixed v,v' and v in M? , we define a unit Jacobian change of variables 

v'l = v'liv[;v,v' ,v) = v[ + ^^^^^{{v - v) ■ ^^^^^}. 
In this change of variables v" € -Ecu" */ cind only if v[ € Eyy' . 
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(40) 



Proof. Assume (39) and (40 1. Clearly 



9"i'(fi) 

dvi 



I. We can check following equality : 



« - V) ■ iv" -v)^ (./ - V + ^^^{{v - V) ■ ^^^}) ■ W - V) 

V — V 

= (vi — v) ■ iv' — v) + \v' — v\{(v — v) ■ -} = (v'-, — v) ■ (v' — v) + (v ~ v) ■ iv' — v) + iv — v) ■ (v' — v) 

\v' — v\ 

= {v'l — v) ■ {v' — v). 

By definition, v'l G E^yi is equivalent to {v'i~v)- ,y,Zy\ = 0- Then, from the above equality, we conclude {v'{—v)- iy,,Z^i = 

which is equivalent to v'{ € Eyy". ■ 

Under the first change of variables ( 39 1 , we can rewrite the first term of ( 38 1 as 

2 



-^—^\wiv")hii,x,v") - w{v')hit,x,v')} f wiv'l)hit,x,v'l)B(2v - v" - v'l, ,""„ ''lA dv'Idv' . 
w — -y'l^ L J J^ ,, \v" — v{\ 



(C) 



(V) 



Under the second change of variables (40), we can rewrite the second term of (38) as 

w{v')h(t,x, v') 



\v - w'l^ 



{£) 



(•^) 



We will estimate (41) and (42) separately in following two sections. 



(41) 



{w{v'l)h{t, X, v'l)B(2v - v" - v'l, fl ""] ) - w{v',)h{t, X, v',)B{2v ~ v' - v',, /^^A } dv'.dv' . (42) 

L ^ \v" — v{\ \v' — v\] } 



3.2 Estimate of ( |41[ ) 

We divide into several cases : 

Case 1 : jwl > A^. From Lemma [6J for A^ > we can estimate 



Q+{wh,wh){t,x,v)l\y\yi^ < C'||/i||^l|^|>Af / _^ w{v') f _ ^ 



,,/|2-7 



dv 



/ 1 1 \ C 

- ^Il^ll- [jiTW ^ (1 + \v\)'-V ^'"'^'^ - N 



Hence we have 



(|4l|l|d>Ar< - 



A^ 



Case 2 : |w| < A^ and \v'\ > 2N, or \v\ < N and \v'\ < ^. Also assume < (5 << 1 



2 X 1|,|<^ / (C) / {V) dv'Idv' 

'{|ti'|>27V or |D'|<i} Je^^„ 



< Cllv 



\v\<N 



< c 



\v'\>2N 



1 



1 



1 



iV2 iV2-7 



+ o{ 



M 



v'|2-7 



2 

oo 7 



k,'|2 1,2 

— I" I , I 2— 

e » dv e ■i 



C 



\v'\<-, 



L. \V 



'|2 |„'|2-7 



(43) 



e <^ dv e * 



oiw) 



(44) 
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where we have used w{v') < e r~ and w{v") < e s~e~ and Lemmapl 
Case 3 : \v\ < N and jj < \v'\ < 2N. 



2 X l|t,|<Ar / 

< c\\h\u 



(C) / (P) dv'{dv' 

\v'\<2N Je^^ii 

1 1 



-ti<\^'\<-^N 



L|u|<A' 



,/|2-7 



|w(i; )h{t,x,v ) — w(v )h{t,x,v )\dv 



(45) 



Since ( i^_^,|2 + |^_^/|2-t ) is integrable we can choose a smooth function z{v,v') with compact support such that 



C\\h\\] 







r 




^ 1 1 A . 






sup / 

|-u|<Af J|t)'|<27V 




,«'J 




V|w-i;'|2 |w-t;'|2-7y -^^^ 


)und ( 


45 


) by two parts 




-/ 




/ 


1 1 \ / ,x 




v'\< 


'^\v\<N 
:2N 


\\v 


-«'P ' |.-.'P-.J ^("'") 


e 



d^;' < 



N' 



-dv' 



(46) 



(47) 



+ C sup |z(u,u')| X ||/i||l=o / l|„|<Ar|w)(w"(f'))/i(i,S,w"(w'))-w(w')^(i,a;,'y')M«'- (48) 

\v\<N,\v'\<2N J^<\v'\<2N 



From (461, it is easy to control the first term 



IHI< §11^11^ 



(49) 



Now we are going to estimate the second term (|48]). Applying Lemma [5 to w{v')h{t,x,v'), we can choose S 
5{wh,n,e,N^,x,jj,2N)>0andaiiopensetUxC{jj<\v\<2N}with\Ua: 



< -^ such that 



\rB{v"iv'))hii,x,v"{v'))-wiv')Ht,x,v')\ < — , 



for v' E {v G W : jj < \v\ < N}\Ux and \{t,x,v) — {t,x,v)\ < S. Therefore we can split the second part (48) as 
integration over Ux and U^ and control as 



C sup |z(t-,t-')|x||/i||^xm3([/,) + C||/i| 

|D|<iV,|i>'|<2Af 



< C sup |z(-i;,-i;')| x 

|ti|<iV,|u'|<2Af 



L^ + CIM^N^^^ 



{T^<|i''l<2A'}n(7- 
3 



\wiv"{v'))h{i, X, v"{v')) ~ w{v')h{t, X, v')\dv' 



(50) 



In summary, combinig (43 1, (44), (49) and (50), we have established 



(41l<C||/i|| 



I -'^ -'W |t)|<Ar,|t;'|<2Af ^^ * 



Choosing sufficiently large N.M > and A''* > then 



(41) < 



C\\h\\^N' 



N^ 



(51) 



3.3 Estimate of (42) 



The estimate of (42) is much more delicate. The reason is that we cannot expect /g (J^) dv'i in (42) is small for all 

w' S K''. We know that h{t,x,v'i) may not be continuous on v'l G &x- Even &x is radial symmetric and has a small 
measure by Lemma pi a bad situation, the intersection of 0^; and Ey^i could have large (even infinite) 2-dimensional 
Lebesque measure, can happen. However we can show that such bad situations only happen for very rare v"s in M.^. 
Using the integration over v' gM.^, we are able to control (42) small. 



Recall {£) and (T) in (42). We divide into several cases : 

Case 1 : |t;| > iV. Follow exactly same proof of Case 1 of the previous subsection, we conclude 



(|42|1|„|>^^ < - 



N 



(52) 
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Case 2 : |w| < A^ and \v'\ > 2N. We go back to original formula, the second term of (38), and use Lemmal6]to estimate 

J\v'\>2N JE„„, J\v'\>2N \V-V\ \I\ 



'|'u'|>2A' •'E^^, J\v'\ 

Case 3 : \v\ < N, \v'\ < 2N, and \v[\ < ^ or \v[\ > N. In the case of \v[\ < j^, we have 
2 X 1^„^^^, 

N 



Ar2-7 



(53) 



i-\v\<N / {£) / (-7^) dv[dv' 

J\v'\>2N J{\v[\<jf}nE^,^,, 



7V^ 



7V2-7' 



In the case of jt^l > iV we have 



2 X l|,|<Ar / {£) / (J-) dv'idw' 

'|i,'|>2Ar J{\v[\<-L}nE^^, 



< 2 



"'^"'^ .d.' 



i,'|2' 



<C||/i|lLe-^ 



{|ti;|>iV}nE„, 



^-^e^{AN+ — + Sy + e-^-^{AN + —y \ dv[ 



(55) 



Case 4 : |u| < N, \v'\ < 2N, and jf < \v[\ < N. In order to remove the unboundedness of |^_"'^,|2 in (42), we choose 
positive smooth function Z{v,v') with compact support such that 



sup 

|-u|<Af J|t)'|<27V 



\v - u'|2 



Ziv^v'] 



dv' < 



m 



Splitting 2 X l|i,|<jv L/|<2Ar(^) /i<l„'l<Ar(-^) dv[dv' into two parts 



<\v{\<I 



2 X l|„|<Ar / w{v')\h{t,x,v)\ 

l\v'\<2N 



-^ -Z{v,v) 



E^,.,n{jf<\v[\<N} 



(56) 



{T) dv.dv' KCWhWl^-^N-'+Mbl) 



ATio 



C 



\v'\<2N I 



sup |Z(w,w')l / {^) dv[dv' 

v\<N,\v'\<2N "'-E„„/n{i<|D;|<iV} 



(58) 



where we used ( 56 ) for the first line. From now we will focus to estimate ( 58 ) . 

Case 5 : |u| < N, \v'\ < 27V, ^ < \v[\ < N and \2v - v' - v[\ < jj ot \v' - v[\ < ^. This region included the part 

where the collision kernel B{-,-) has a singular behavior. 

CI \\h\\oc sup \Z{v,v')\ {J") 'i-{\i2v-v')~v[\<^ or \v'^v\\<i,}iv',v[) dv[dv' 

\v\<N,\v'\<2N •'E^^,n{jf<\v{\<N} ' " ' '■' 



'|u'|<2Af \i 

< C sup |Z(t),w')| X 11^112 

|i>|<Ar,|D'|<2Ar 

< C sup \Z{v,v')\x 



dv' 



|'u'|<2Af 



N^ 
'M2- 



'e '"4' f dz;i{l{|(2„_,,)_,,|<^}(^'l) + l{|,,„„-|<^}(i;l)}xiV^ 

(59) 



|D|<Af,|i;'|<2Af 

Case 6 : \v\ < N, \v'\ <2N,j.< \v[\ < N and \2v - v' - v[\ > ^ and \v' - wil > ^^ and < (5 < 3^. We estimate 



2 X l|i,|<jv / dv''w{v')h{t,x,v')Z{v,v') / {u;(u")/i(i,a;,'y")-S(2w - w" - u", -r^; ^) 

'|D'|<2Af "'S„„,n{i<|tii|<Af} 1^' - «1 I 



-«)(«;)/i(t,a;, v[)B{2v - v' - v[, y-^^)}l{^2v-v'-vi\>i,}h\v'-v[\>^}dv[. 



(60) 



We need this step because of the singular behavior of 



B{ui,U2) = \uipqQ{- — - ■ - — -) = \ui\'^{qoo^){ui,U2) 

\ui\ \U2\ 
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where -5 : M'^ x M'^ — > M with ^{ui,U2) — -njn" ' nfi- The function 5'(wi,U2) is not continuous at (ui,U2) = (0,0) and 
continuous away from (0, 0), i.e. the restriction of 5^ on a compact set, 

dM,N : {^ < l^il < 6iV} X {^ < \U2\ <AN}^R 
is uniformly continuous. From \2v — v' — v'i\ > jj and |f — u| < 5 < j^j^ we have lower bound of 



\2v - v" - v'l\ > 



lOAf 
V' — V 



V — V 



\2v ~ v' - v[\ ^ \v - V ^ j-,—^{{v - V) ■ j-,—^}\ 



> 



2M' 



Similarly from \v' — v'i\ > -^ and |w — wj < 5 < j^ we have a lower bound of 



Therefore for any e > 0, we can choose (5 > so that 



> 



2M' 



B{2v - v" - vl ^7^^) - B{2v - v' - v[, ^^^-\) 
\v" — v'-[\ \v' — V-y\ 

\2v - v" - <r((7o o :S){2v - v" - v'i, v" - v'i) - \2v ~ v' ^ v[r{qo o 5)(2« ^ v' - v[,v' - v[)\ < ^, (61) 



for \2v -v'- v'i\ > jj and \v' - v{\ > jj and < J < j4 



lOM • 



Now we split ( 60 ) by two parts 

2 X l|„l<Ar / dv'.. / w{v'()h{t,x,v'l)\B{2v~v" -v",-— jjr) - B{2v - v' - v[, — ]-) 

J\v'\<2N JE^^,,n{j,<\v[\<N} { \v"-v'{\ \v'-v[\ 

+ 2xl|,|<^/ dv'.. f {w{v'l)hii,x,v'l)-wiv[)hit,x,v[)}Bi2v-v'^v[,^^j^^). (62) 



Using (61 1, the continuity of B{-, •) away from (0, 0), the first fine above is bounded by 



Csup\Z{v,v')\ X 



'N, 



(63) 



In the remainder of this section we will focus on ( 62 ) 



Estimate of (62) 



p2\ <CN^\\h\\^ sup \Z{v,v')\ I w{v') I \w{v'l)h{t,x,v'l) - wiv[)h{t,x,v[)\dv[dv' , (64) 



where we used sup|^|<^ |,„,|<2Ar.|t,'|<Ar B(2v — v' — v[, Z,_^Ji ) < oo. Recall our choice of v" and v'{ from (39) and (40) to 
have 



\v'l-v[\< 



T-, — r{(^ - '^) • r^ — f} <\v-v\<5. 

\v' — v\ \v' — v\ 



We 



will use the foUowin strategy : separate /g ni ^<\v' \<n\ ■■■dv'i into two parts 



/ ...dv[ + I ...dv[. 

The first part is the integration over Ux, a neighborhood of iSx that contains possible discontinuity of h. Moreover we 
expect the measure of the neighborhood Ux is small so we can control the first term. For the second term, we will use 
the continuity of the integrand wh. However if w = then 0x could be a large measure set in Ey^i n {jj^ < |wi| < N}. 
For example if ©i: n S^ = {w G §^ : "3 = 0} then (&x is xy— plane and Eqs^ is also xy— plane. Therefore we have to 
divide two cases v ^ and v ~ and study separately. 
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Case of V 7^ 

In the case oi v ^ 0, assume g < \v\'^/2 for sufficiently smaU g > 0. We will divide the velocity space M.^ into 



m ^ Iv' eR^ : \v\ - A < v' ■ ^ < \v\ 



and ^^ ^ Iv' e 



1 V 






- \v\ 


\v\ 





> 



(64), ♦ into 



The important property of *B is that if w G S"^ then E^yi does not contain zero. We can split the integration part of 

X 



w{v) / \w{vi)h{t,x,Vi) — w{v[)h{t,x,v[)\dv[dv' 

'eS2jvn<8 JE^^'n{jf<\v[\<N} 



w{v') / \w{v'l)h{t,x,v'l) - w{v[)h{t,x,v[)\dv[dv' . 

v'€B2n\'B "'-E„„,n{i<|t;;|<Af} 



(65) 
(66) 



Notice that 25 H Sow has a small measure 



m3(S n B2n) < 2TT{2Nf x 2 j^ < 27r(27V)2 x 2-^ < 2V2TT{2Nf^. 



Therefore we have 



\[(Sh)\<CN^\\h\\L^^g. 



(67) 



Now we are going to estimate ( 66 1 . Here we use a property of QS'^ : for v' e *B'^ we have 

\v' ■v-\v\^\ 



dist(0,-Ei,i,') 



> 



> 



> 



\v' -v\ \v' -v\ 2iV+|w| - 37V' 



where we also have used |w'| < 2N and |i;| < TV. From Lemma [5] we use U^^ an open radial symmetric subset of 
{h ^ l"!^].! !i N} with a small measure and wh is uniformly continuous on C/^, to split (66 1 into 



w{v') / \w{v'l)h{t,x,v'l) - w{v[)h{t,x,v[)\dv[dv' 

v'<£B2n\'B "'s„„'n{i<|i.;|<Af}n(7^ 

w{v') / \w{v'^)h{t,x,v'l) - iu{v[)h{t,x,v[)\dv[dv' . 



(68) 
(69) 



For the last line, we use Lemmalslto know estimate \w{v'{)h{t,x,v'{) — w{v[)h{t,x,v[)\ < -^, for v[ € Eyy' H {-^ < 



\v'i\ < N}\U^ and [< - v[\ < \v -=^\ < S. Therefore 



(69)| <CN^ 



N, 



(70) 



In order to show that (681 is small, we introduce following projection : 

Lemma 9 Assume < g < ^. Let Ey^i = {v[ G M"^ : {vi — v) ■ [v' — w) = 0} . We define a projectic 



§ — > E„„/ 



M e S^ ^ 



V ■ (v' — v) 
u ■ (v' — v) 



u € Ey 



For v' G {v' G M'' : |v'| < 2A^}\*B, define the restricted projection 

P' = P\p-iiE^^,,n{i/N<\v[\<N}) ■■ P'HEyy, n {1/iV <\v[\< N}) -^ Eyy, n {1/iV <\v[\< N}. 
Then for v' G i?2Ar\23 the Jacobian ofV is bounded : 

3N'^ 



Jac(P') = 



where 9 is defined by cos 6 — u ■ |",_^| . 



9m 



w — v\ 



sec^6ltan6' < 



Proof. Without loss of generality, we may assume -rjrz^i — (0,0, 1)"^. Using spherical coordinate. 






v' -V I sin B cos 6 \ , I tan Q cos ( 



u ■ (v' — v) 



\v' — v\ 



cos( 



sm t/ sm ( 
cos6' 



V ■ 



V — V 



\v' — v\ 



tan 6 sin ( 
1 
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Figure 4: Projection Map 



and a Jacobian matrix of P' , 

Therefore a Jacobian of P' is 

Jac{¥') 

Notice that 

|sec6'| = — 



d¥' 



V ~ V sec" u cos ( 

'" ■ — ■ I „2 , 



tan 9 sin ( 



d{9,4>) \v' ~ v\ \ sec 61 sin ^ tan 6* cos < 



dP' 



d{e,.p) 



sec^ 6*1 tan 6*1 < dist(0, ^1,^,0^1 sec 6*1^ 





1 








v' 


— V 




'(/, • 










\v' 


—v\ 





V - {v' ~ v) 
u • {v' — v) 



V • 


v' 


— V 


\v' 


-v\ 



\r{u)\ 



dist(0,£'„„') 



Because V'{u) e {^ < \v[\ < N} and dist(0,^„t,') > gf^ we have 

|P'(u)|3 3N^ 



|dist(0,£;w)| " g 



Assume we choose m2{Ux n S^) < j^jji- By definition we know that V'{Ux n S^) = E^^, ^ {j; < Wi\ < N} r\Ux and 



the 2-dimension Lebesque measure of E^yi H {-^ < \v'i\ < N} D U^ is bounded by 



SiV" 



m2(^wn{- < \v[\ <7V}n(7,)=m2(P'(C/,n§"))< Jac(P')x |c/,ns'| < X 



-N 



3N' 



g N,N^ gN, 



e. 



Therefore we have an upper bound of ( 68 1 



(68 1| <CN^e 



where C = /^a w(v')dv' . In case oi v y^ 0, from (67), (70) and (71 1, we have 



(621 <CN^\\h\\^ sup \Z{v,v')\ x i < CN 



4|| J,l|2 



^ sup \Z{v,v')\{N'^g+{l + )— }. 

\v\<N,\v'\<2N Q ^^* 



(71) 



(72) 



Case of V = 

In this case, we do not have a upper bound of the Jacobian of P'. Instead we will use the structure of i5x of Lemma [1] 
crucially. In case of i; = 0, we split (64 1 



w{v') / \w{v'{)h{t,x,v[) - w{v[)h(t,x,v[)\dv[dv' 

\v'\<2N JEa^,n{-L<\v[\<N} 

w{v') I \w{v'{)h{i,x,v[) - w{v[)h{t,x,v[)\ X luAv'i)dv[dv' 

\v'\<2N JEg^,n{jj<\v[\<N} 



(73) 



w{v') / \w{v'l)h(t,x,v[)-w{v[)h{t,x,v[)\xlE ^^i^^^,,^^y^^jj{v[)dv[dv' (74) 



For v', we use a spherical polar coordinates (r',6'', (/)') so that 

v' = (r'sin0'cos0',r'sin6'' sine/)', r'cos^'). 



(75) 
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By definition, Eq^i is a plane containing the origin and normal to v' . We know that Eqv' is generated by two unit vectors 



Eqv' 



cos f COS ( 

cos 9' sin ( 
-sin 61' 



— sin( 
coscj) 




We will use a polar coordinate {r'i,9'i) for v'l e -Eou', i-e. 



{v[)i \ I cos f cos 0' —sa\(p' smf'cost 

K)2 \{r[,e[-e',(l)') = r[\ cos 6*' sin f cosf sin6''sin0' 

M)3 







cos( 



COSW^ 

sin e'^ 





Direct computation gives det 



[dir[.e[.e')) - 



sm y cos 
Therefore we have following identity 



smu smf 



cost 



p poo pZTT PTT 

■■■dv[^ / / ■■■ {r[f cos e[d9'd0[dr[. 

JR^ Jo Jo Jo 



Recall the standard 3-dimensional polar coordinates and 2-dimensional polar coordinates : 

(.2N r2TT n-n 



■■■ dv' ^ 



|i;'|<2Af 



Eo^'n{Jr<\vi\<N} 



/ / / ■■■ {r'fsme'de'd(t)'dr', 
Jo Jo Jo 



N I.21X 



■"'-/J 



• r[d0[dr[, 



and use above identities to control ( 73 ) by 

r2N /.27r 



N 



Kl! 



2 /•2-n- 



dr'{r'yiD{r') / d0' / d6''sin6l' / dr[r[e ^ / de[luAv'i{r[,S'i;0' ,(/)')) 

Jo Jo Ji Jo 



(76) 



cos a cos cos tli — sm (p sm tla[ — cos cos </> sm t^j^ — sm cp cos t^j^ sm cos 
(rj^)^ cos 6*5^ det I cos 0' sin 0' cos 6*']^ + cos 0' sin 6*5^ — cos 6*' sin 0' sin 615^ + cos </>' cos 6*^ smd'sincj)' I = (r'j^)^ cos^'j^. 



(77) 



(78) 



We focus on the inner integration (^ and divide it into 






d^'sin^' I dr[r[e-^ / d0;ic/jz^;)l,ie(f-,.f +p)u(^-e.^+.) 



/■TT pN , I ,2 /'2'7r 

/ de'sinO' I dr[r[e-^^ de^luAv'i)^', 



;e[o.f-e]u[f+e,Tf-e]u[4f+e:27r]- 



1 JV^ 



and r']^ e [j^,N] to have 



(79) 
(80) 



Easily (|79|) < 2g{e-^^ -e-"^) < ig. For (|80|), we use 1 < ^^^^ and 1 < Nr[one[ e [0,^-g]U['^+Q,^-g](J[^+g,27r] 



TT pN 



2tt 



{80)<g-'N dd' dr[{r[f dd[ cos 9[lujv[{r[, e[; 0' ,(j>')) ^ g'^N x m^iU^n {-< \v[\ < N}), (81) 



N 



where we used ( 77 1 . To sum we have 



(73) < (781 < C\\h\\^ Ug+g-^N x ^| 



(82) 



On the other hand for ( 74 ) we can use Lemma Is] to have 



(741 <C 



N, 



(83) 
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From (182j) and (|83|) we have 

(62|<C7V2||/i||^sup|Z(i;,i;')| x4 



CN' 



,sup|Z(i;,i;')|x{(73l + (74)} 



< CN^\\h\\^snp\Z{v,v')\l^ + \\h\\^[4g + g-'Nx^] 



(84) 



To summarize, from (|52|), (|53|, (|54|), (|55|), (|57|), (|59|), (|63|), (|72| and (|84|), we have established 

sup 



1 



JV^ 



|42|)<C||/i||^{-+e-T^} + q|/i||^ sup \Ziv,v')\—+C\\h\\i^ 

^^ |i>|<Ar>'|<2Af ^^^ \v\<N,\v'\<2N 



N, 



■c\\h\u 



sup 

\v\<N,\v'\<2N 



Z{v,v')\{N^^ + AN^g) 

<2N 

\Z{v, v')\ U^ + \\h\\^(l + ^+N^ + N^g^ I . (85) 



We choose N, M, N^, > sufhciently large and g > sufficiently small so that we can control ( 42 1 < | . Combining with 
the result of previous subsection (51 1, we conclude (37 1 and and prove Theorem El 



3.4 Continuity of Collision Operators Kf and T(f,f) 

The following is a consequence of Theorem |4] 

Corollary 5 Assume f{t,x,v) is continuous on [0,r] x (f2 x IR'^)\0 and 

w{v)f{t,x,v) - {l + p'\v\yf{t,x,v) e L°-{[0,T] X (n X M^)). 
Then K f{t,x,v) andT^{f,f){t,x,v) are continuous in [0,T] x fJ x M"^ and 



sup 1^ ^{v)w{v)K{f)\ < 



[0,T]xnxI 



sup W-\v)w{v)T+{fJ)\< 

[0,T]xOxR3 



Proof. The above boundedness is a direct consequence of ( 33 1 and ( 34 ) . Thanks to Theorem l4J we already established 
the continuity of r+ . Therefore we only need to show the continuity of 

——Q^{^/Jif,li) — e * / / B{y — u,uj)f{t,x,u)e ^ dcjdu. 

Choose {i,x,v) ^ {t,x,v) so that \{t,x,v) ~ {t,x,v)\ < S. We will estimate 

-—Q^{^f,fi){i,x,v) —Q^{^f,^i){t,x,v) 



1 



S2 



{B{v — u,uj)f(t,x,v) ~ B{v ~ u,Lu)f{t,x, u)}djjjdu 



1 



. . B{v — u,uj)e * f(t,x,v)diodu — / / B{v — u',uj)e * f{t^x,u' — {v — v))dujdu' 

VM JK3 Js2 VM JR3 Js2 



< —=i I I \B{v-u,uj)\\e 4 _e ^ \w ^ {u ~ {v - v))\\w f 

Iff _M^ 

— / / \B{v — u,uj)\e " \f{t,x,u) — f{t,x,u—{v — v))\diodu, 



-jdujdu 



(86) 
(87) 



where we used a change of variables u' = u + {v ~ v) for the underlined term. Using the Taylor's expansion we control 



e 4 = e 



1 l".|2 1 a^ |u|2 

+ x|w*|e ^1^ - ^1 < -(|u| +^)eTe~T^ x |w - wl, 



where w* = s*{u — (v — w)} + (1 — s^,)u for some s.^ G (0, 1) and |w — f | < S. Therefore we control 

.1 



|(|86l)| < e' i I \v ~ up ^{\u\ + S)e'i e ^ du x sup qo(- r ■ uj)duj\v - v\ x \\wf\\ao 

2 v,u Js^ ^\v-u\ 



< C{l + \v\re—\\wf\U 
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where we have used the the angular cutoff assumption dsl). Now we estimate (871 with following steps : 

„|ii _jv2 _M^ 

Case 1 : |u| > N. Since e ^ < e » e s^we estimate 

/ / :i,dujdu <Ce-^ j e-^-^\u-v['dux\\wf\\^<Ce-^v{v)\\wf\\^. (89) 

Case 2 : |u| < A^. A function / is continuous on [0,T] x (fj x i3(O;Af))\0. By the Lemma p] we can choose 
U^. C B{0;N) with |f/^| < § with \fit,x,u) - f{t,x,u- {v - v))\ < ■§ for \{t,x,u) - {t,x,u- {v - v))\ < S with 
u e B{0;N)\Ux- Therefore /|„|<jv /s2 ^ dtodu is bounded by 



/ / Jltdujdu+ / / Jit dujdu < C—i'{v)\\wj 

Ju&B(0;N)nU^ Js^ Ju&B(0:N)\U^ Js^ ^ 



(90) 



From (881, ( 89 ) and ( 90 ) , we summarize 



^|g_(^/,/i)(f,x,z;) - Q^i^f,f,)it,x,v)\ < io{6) + e"* + ^)'^\\wf\l 
which is less than e for sufficiently large N and sufficiently small d. ■ 



In following sections, we will prove Theorem 1, Theorem 2 and Theorem 3, for each boundary conditions. In order to 
write theorems in the unified way [13] for all boundary condition cases, we use the weight function w(u) = {1 +^^11)^}^ 



in ( 12 ) and define 

h = w[v) X 



In terms of /i, the Boltzmann equation M is equivalent to 

{dt+v-Vx + v-K^}h^wT{-,-), (91) 

w w 

where K^h = wK{ — ) with boundary conditions : 

1. In-flow boundary condition : h(t,x,v) = w{v)g(t,x,v) for (x,v)^^^. (92) 

2. Diffuse boundary condition : 

h{t,x,v) = w{v)y^fi{v) / h{t,x,v')————^=Cf_i,fi{v'){nx-v'}dv' for (a;,u)e7-, (93) 

Jv'-n{x)>Q W[v')y'fl{v') 

with a normalized constant c^ in (rol). 

3. Bounce-back boundary condition : h{t,x,v) = h{t,x,—v) for {x,v) E 7_. (94) 

4 In-Flow Boundary Condition 



In this section, we consider the linearized Boltzmann equation (91) with the in-flow boundary condition (92 1. First we 
will show the formation of discontinuity using a pointwise estimate of the Boltzmann solution 13J. Then we use the 
continuity of collision operators, Theoreml4] to show a continuity of solution on the continuity set £ and the propagation 
of discontinuity on the discontinuity set 2). 

4.1 Formation of Discontinuity 

We prove Part 1 of Theorem [I] Without loss of generality we may assume xq — (0,0,0) and vq — (1,0,0) and 
(xojWq) G 7o- Locally the boundary is a graph, i.e. ft Ci B{0;S) = {(xi, a;2, 2:3) € B{0;S) : 2:3 > ^{xi,X2)}. The 
condition (a;o,t;o) G 7o implies tb(xo,wo) 7^ and tb{xo,—vo) ^ which means $(^,0) < for ^ e (— (5, (5)\{0}. (See 
Figure 3) 
For simplicity we assume a zero boundary datum, i.e. g = 0. From Theorem 1 of [13 , we have a global solution of the 



linearized Boltzmann equation (91 1 with zero in-flow boundary condition, satisfying 



sup \\h{t)\\^<C'e-'^'\\ho\\c 
te[o,oo) 
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for some A > 0. In the proof wc do not use the decay esthxiate but just boundedness 

sup ||/l(t)||oo<C"||/lo||oo. 

te[o,oo) 



Recall the constants Ck and Cr from (33 1 and (34). Choose ig G (0,min{ 



S tb{xo,-vo) 



}) sufficiently small so that 



1 < (e-''(i)*° - ioCkC" - (1 - e^''(i)*°)Cr(C')') , 
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(95) 



(96) 



where ^{l) = j^(uo) for any vq E M.^ with |wo| = 1. This choice is possible because the right hand side of (96) is a 
continuous function of tg G M and it has a value 1 when to = 0. Furthermore assume a condition for our initial datum 
ho : there is sufhciently small 6' — 6'{n,to) > such that B{{—to,0,0);S') C ft and 



/io(a;o,i'o) = ||/io||oo>0 for {x,v) e B{{-to,0,0);5') x B{{1,0,0);6'). 



(97) 



We claim the Boltzmann solution h with such an initial datum ho and zero in-flow boundary condition is not continuous 
at {to, xo, vo) = {to, (0, 0, 0), (1, 0, 0)). We will use a contradiction argument : Suppose 



[H'to)]xo 



0. 



(98) 



Choose sequences of points {x'^,v'^) = ((0,0, i), (1,0,0)) and {Xn,Vn) = ((^,0, 4>(i, 0)), (1,0, ^)). Because of our choice, 
for sufficiently large n € N, the characteristics [X{0; to, sq, wo), V(0; to, Xq, vo)] is near to {{—to, 0, 0), (1, 0, 0)), i.e. 

{x'„ - tov'„, v'J = {{-to, 0, -), (1, 0, 0)) e B{{-to, 0, 0);d') x S((l, 0, 0);d'). 

n 

Hence the Boltzmann solution at {to,x'^,v'j^) is 

h{to,x'^,v'J = ho{x'^-tov'^,v',,)e-'^(<^''>+ f\-i<)i*o-r) Ij^^^^^^y^h^h^ 

Jo I WW) 

= ||/io||ooe-''(<)*«+ [ "e-'^^^^'"-^^ { K^h + wr{-, -)] {r,x - v:^{to - T),v'JdT. 
Jo I WW) 



(99) 



Combining h{to,Xn,Vn) — w{vn)g{to,Xn,Vn) — with ( pJ8| ), we conclude 

/i(tg, x'^, v'J ^ asn^ 0. 



On the other hand, using (95) we can estimate 

liminf |/i(to,a;J^,w^)| = limmi\h{to,x'„,v'J - h{to,Xn,Vn)\ 



> liminf |||/io||ooe~''(''")*" 



Ci^C'WhoWoodT- 



,^{v'Je-''(<^^'-^^Cr{Cr\\ho\\lodr\ 



> ll/^olUe-'-W*'' - ioCkC'11/.olU - (1 - e-''^'^'o)Cr{C'f\\ho\ 

\\h 



2 

oo 



= \\h 



oiloo (e-'^d)*" - toCkC - (1 - e-''(i)*«)Cr(C")') > 



Oiloo 



^0, 



which is contradiction to (99). 



4.2 Continuity away from 1) 

We aim to prove Part 1 of Theorem [3] in this section. First we recall Lemma 12 of [T3j, the representation for solution 
operator G{t, 0) for the homogeneous transport equation with in-flow boundary condition : 

Lemma 10 USf Let ho G L°° and wg G L°°. Let {G{t,0)ho} be the solution to the transport equation 

{dt + v-V,}G{t,0)ho^O, G{0,0)ho^ho, {G{t,0)ho}^_^wg. 
For {x,v) ^ 7on7_, 

{G{t, O)ho}{t, X, v) = lt-ti,<oho{x - tv, v) + lt-ti,>o{wg}{t -tb,x~ t^v, v). 
Next we prove a generalized version of Lemma 13 in jl3j. 
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Lemma 11 (Continuity away from 2) : Transport Equation) Let fl be an open subset o/M"^ with a smooth bound- 
ary dH, and an initial datum hQ{x, v) be continuous in fl xM.^ U {7_ U 7+ U 7o~} o,nd a boundary datum g be continuous 
in [0, T] X {7_ U 7q~}. Also assume q{t, x, v) and (j>{t, x, v) are continuous in the interior of [0, T] x fi x K"^ and satisfy 
suprQ T]xf2xK3 9(^1 2;, u) < 00 and supfQ rjxn k'('j 'j "J^) < 00 for all f G M^. Let h{t, x, v) be the solution of 

{dt +v -Vx + <t>}h = q , h{0,x,v) = ho, h\^_=wg. 

Assume the compatibility condition on j^ U 7q~, 

ho^Xjv) — w{v)g{0,x,v). (100) 

Then the Boltzmann solution h(t, x, v) is continuous on the continuity set C Furthermore, if the boundary dQ does not 
include a line segment (Definition]^ then h{t,x,v) is continuous on a complementary set of the discontinuity set, i.e. 
{[0,T] xClx M3}\D. 

Proof. Continuity on {{0} x 0, x W^} U {(0, 00) x [7- U 7o ]} is obvious from the assumption. Fix (t, x, v) e €. Notice 
that 

j-Jh{s,Xis),Vis))e-i:'f'^^'^(^^'''(^^^^^}-q{s,Xis),V{s))e'f'^<^^'^^^^^^ (101) 

along the characteristics X{s;t,x,v) — x — v{t — s),V{s;t,x,v) — v until the characteristics hits on the boundary. 
Choose (i, x, v) ^ (i, x, v) and use a change of variables s = s — {i — t) with s e [t — i,t] to have 

\-{h{s + (f- t),X{s), V{s))e- f^ Hr+ii~t),Xir),Vir))dry 

Ids 

-9(5+ (t~t),^(s),t^(s-))e-/^^(-+(*-*)'^(-)'^(-))'*-}l[_(,_,)+,„,,{o,*-t.(s,.)},t](s) - 0, (102) 

where X{s) ^ X{s+ {t ~ t);t,x,v) and V{s) = V{s + {t - t);t,x,v). 

By the definition £, we can separate two cases : t < t)3{x,v) , {x-h{x,v),v) G 7_ U7q~. 



Case oit — tb{x,v) < From the assumption t — tb{x,v) < 0, we know that (101) holds for < s < t. Now we 
choose {i,x,v) near {t,x,v) so that t — t-h{x,v) < 0, and X{s) — X{s + {t — t);t,x,v) is in the interior of il for all 
s £ [t — i, t] . Taking the integration over [min{0, t — t},t] of ( |101[ ) — ( 102 ) to have 

h{t, X, V) - hit, X, V) = /lo(^(0), l/(0))e- J'o 4>{r,X(r)y(r))dr _ ^^(^(^ _ ^ ^y ^^ _ ^-^^g- //_ ^ 0(r+(t-t),X(r),y(r))dr 

+ /* |l[„.ax{0,t-t.(.,.)},t] (S)g(s, X{S), V{s))e- ^l ^ir,nr)y(r))dr 

Jmin{0,t-t} ^ 

-l[t-f+max{0,t-t.(2,C)},*](5)9(s+(t-t),X(s),t/(s))e-/>(-+(*-*)'^M'^M)'^-}ds. 



Since /iq and </) is continuous, it is easy to see that the first line above goes to zero when {t,x,v) — > (t,x,v). For the 
remainder we separate cases : i — i > and t — i<0. Ifi — t>0 the remainder is bounded by 

|<j(s)er *M- - q(s + (ri - t))e--''X"+(*-*)| + |t - t\ sup ||(7(s)|Ue*™Po<=<' 11^(^)11-, 

0<s<t 

where the first term is small using continuity of q and 0, and the second term is small as ij,,x,v) — >■ {t,x,v). The case 
i — i < is similar. 

Case of [xy,[x,v),v) e 7_ U 7q^ We only have to consider cases of t > tTa{x,v) and t = t\,{x,v). By definition 

(a:b(a;, v), w) S 7_ U 7q". From Lemma^l we know that t\y(x, v) is a continuous function when [x\,(x, v), w) ^ 7_ U 7o~- 
In the case of t > t];y(x,v), for {t,x,v) '^ (t,x,v), we have t > tiy{x,v). Taking the integration over [min{0,i — i},t] of 
( [101] ) - ( [T02| ) to have 

hit, X, v) ~ hit, X, v) = wgit - tUx, v),Xity,ix, v)), Vit^ix, v)))e- ^-*.(==-') ^(-'^M.^M)'^- 

-«;5(t-ib(x,f)),X(tb(S,f;)),F(ib(S,i;)))e--^''-'.(^.^)^(^+('-*)'^(^)-^(^»^^ 

9(s,X(s),y(s))e-/>(^^^(^)'^(-))'*^ds 

9(s + (t-t),X(s),\/(s))e-^='*(^+(*-*)'^(^)'^("»'*^ds. 
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Using the continuity of ib and q and 0, it is easy to show that \h(t,x,v) — h{t,x,v)\ — > as (t,x,v) — >■ {t,x,v). In 
the case of t — tb{x,v) we can choose (t,x,v) ~ (i,a;,w) so that tb{x,v) € {t — e,t + e). Taking the integration over 
[minjO, t~t},t] of ([lOl|) - (Il02]) to have 



|/i(f, X, t;) - h{i, x,v)\< wg{t ~ %{x, v),X{t^{x, v)),V{t^{x, v)))e- -^-'.t-^) ^(-■^M^^M)'^- 

g(s, X(S), y(s))e- ^s 4>ir,X(r)^V(r))dr _ ^(^, ^ (^-_ ^)^ ^(^)^ ^(^))g- /^* 0(r+(*-t),X(r),yM)dr 
t — tb{x,v)+e 

+2e sup ||g(s)||ooe*™P''2=^*ll'^(")ll~, 

0<5<t 

where the first three hncs can be small using the compatibility condition and continuity of /iq in rj x M"^ U {7- U7+ U7q^} 
and a continuity of g on [0,r] x {7- U 7q~} and continuity of <j). For the fourth line above, we use the continuity of q 
and (f). 
If the boundary dfl does not include a line segment (Definition^ we have € — {[0,T] x Q x R^}\'D. ■ 

Proof of Part 1 of Theorem |3] 

We will use the following iteration scheme 



ds 



{dt + v-V^ + v}h 



m+l 



K„X 



wFj 



'^m ujn 



W W 



wT^ 



f^m+1 



w w 



/ith h"'^^^\t^Q — ho and h^^^^{t,x,v) — wg{t,x,v) with {t,x,v) e 7_ U 7q . For simplicity we define 



q" = if^/t" + wT^ 



/i" h" 



— wF_ 



h"' h"'+^ 



(103) 



(104) 



Step 1 : We claim 



for alH e N and for any T > where 



/i* is a continuous function in Cy 



e:n{[o,T] x r? x m^}^ 



(105) 



(106) 



(1051 is satisfied for i = 0. Assume (|105 
equation of /i™+^ is 



for all i = 0, 1, 2, ..., TTi. Rewrite wF_ 



w ' w 



wher e the continuity set £ is defined in ( 10 ). We will use mathematical induction to show ( 105 ). We choose h'^ ~ then 

= ly (yMiJr) /i™+^ then the 



{dt + v-V^ + v{v) + V 



w 



771 + 1 



K„h" 



- wFa 



(107) 



From Theorem Hand Corollary p^ we know that v (y^^) and wF+ (^, ■^) is continuous in [0,T] x fi x M^. Apply 
Lemma 11 where (f>{t,x,v) corresponds to v^v) + vi^-s/v^) and q{t,x,v) corresponds to the right hand side of (107 1. 
Then wecheck ( |105[ ) for i = rn + 1 . 

Step 2 : We claim that there exist C > and (5 > such that if C{||/io||oo+supo<<;<oo ll'"^5(5)lloo} < ^ and C||/io||oo < 5 
then there exists T == T{C, (5) > so that 



sup \\h^{s)\\^<C\\K\ 

0<s<T 



(108) 



for all TO e N. Moreover {ft.™}^=o is Cauchy in L°°([0, T] x Q x E^). 

First we will show a boundedness ( 108 1 for all to e N. We use mathematical induction on m. Assume supg<- j<-rp 1 1 ft,™ (s) | |oo < 



C'llftolloo where T > will be determmed later. Integrating (103) along the trajectory, we have 



' ?n+l 



{t,x,v) = lt<:ti,(x,v)e '''^''^*ho{x-tv,v) + lt>t^,(x,v)e ''^''^*''''''''"''w{v)g{t - tb{x,v),Xb{x,v),v) 



m.ax.{t — t}:i{x,v),0} 



"'(^)(t-^){K^h'"' + wT+ 



I ?n L7n 



wT^ 



ft™ ft™+l 



} (s, X — (i — s)w, w)ds 



< \\ho\\oo + sup ||wg(s) 

0<s<t 



iCk sup ||ft™(s) 

0<s<t 



Cr sup ||ft'"(s)|U sup (||ft™(s)||^ + ||ft™+i(s)||^) 

0<s<t 0<s<t 
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and 

sup \\h + (s) oo < z ^ r<fiih II — Tl n f Ml \ i K^o 00+ sup \\wg{s) 

Q<s<t 1 -CrC{||/io||oo +supJ|wg(s)||oo} L o<s<t 

< C < \\ho\\oo + sup ||wg(s)||, 
L o<s<t 
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C-3 



where we choose C > 4 and then {||/io||oo + supo<s<t 11^3(3)1100} < 2cr~c ^^'-^ O^en ^ — 2c c 

The equation of /i™+i - /i™ is 



= 



Newt we will show the sequence {/i™} is Cauchy in L°°([0,T] x fl x 



rn+l L'^n\ 



(/i™+^-;i")|t=o = o, ih"'+'-h 



(109) 



where 



q^ = K^ih"" -ll^ ^) + tor4 



/l'" ^m _ ^T7 



u; w 



-Ji.Fa 



-ior_ ( — , I +ior_ 



w w 



w w 



(110) 



From ( 33 1 and ( 34 ) , we have a bound of g,„ , 

sup ||r(s)||oo<Ck sup ||{/i"-/i™-^}(s)||c 
I 771 1, rn — 1 T 



(111) 



+Cri^(i;){ sup ||{/i'"-/i"-^}(s)|U+ sup ||{/i™+^ - /i"}(.s)|U} x ( sup ||/i™(s)|U + sup ||/i'"+i(s)|U). 

0<s<t 0<s<t 0<s<t 0<s<t 



Integrating ( 109 ) along the trajectory, we have 

nt 



||{/i'"+i-/i'"}(i)|U < / e"'^^^^'-^^\\r{s,x-{t-s)v,v)\Uds 



< Cki sup \\{h"'-h"'-'}{s) 

0<s<t 



+CCr ||/io||oo+ sup \\wg{s) 

\ 0<s<t 

If we choose CCr||/io||oo < \ and CkT < \ then 



sup ||{/i™-/i"-i}(s)|U+ sup \\{h^^'-h^}{s)\\, 

0<s<t 0<s<t 



sup ||{/,™+^-ft™}(s)||^< sup ||{/i™-/i"-^}(s)||, 

0<s<T ^ 0<s<T 



Then we have 



sup ||{/i'"-/l"-n(s)|U < sup ||{ft'"-/i'"-n(s)||oo + ---+ sup ||{/l"+^-/l"}(s)||oo 
0<s<T 0<s<T 0<s<T 

< i^^^^ + --- + i} ^^P \\W'+' - h'^mW^ < ^ sup ||{/^l-/^°}(.) 

2™ " 2" 0<s<T 2" o<s<T 

< 7^C{\\ho\\oo + sup ||«;5(s)||oo}, 

2" 0<s<T 



which means that the sequence {/i™} is Cauchy in L°°([0,T] x f2 x 
Step 3 : From previous steps we obtain that h with lim„_>oo h"' is continuous function on £7. Now we claim that h is con- 
tinuous in £. Notice that T only depends on ||/io||oo andsupo<s<2n ||w(7(s)||oo- Using unform bound of supQ<^^o^ ||/i(s)||oo 
(Theorem 1 of [T3|) we can obtain the continuity for h for all time by repeating [0,T], [T,2T], .... If the boundary dH, 
does not include a line segment (Definition p| then every step is valid with [0, 00) x {Jl x K^}\S) instead of € and 
[0,r] X {nx M3}\J) instead of €t. 



4.3 Propagation of Discontinuity 

Proof of 1 of Theorem [2] 
Proof of ( [18| ) 

In order to show the upper bound of discontinuity jump (18 1, we will show 

\h(f)] , ^ <\h] p-{Tt + -Sij\\'^o\\^){l + \vo\rit-to) 

['l{t)\xo + {t-to)vo,vo ^ ['T'\to.xo,voe ^^ "^^ ) 



(112) 
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when (^ojiio) e 7^ and t G {to, to +tiy{xo, —vq))- Choose two pohits {x',v'), [x" ,v") € {fj x ]R^\©} n B{{x,v);5)\{x,v) 
and compare the representation 

\h{t,x ,v ) — h{t^ X . V )| < lt-tn>t ix' v')hii ~ tbi^ ^ ^ )t ^hi'X T V ), V )e *-*b(^'^) ^ 

+ ^t-to<t^^^x'.^^)h(to■.x '(t~to)v ..v)e "> '*o ^Va-^a. 

„ „ ,, „ „ -I/(u")t|3(^",l>")-/'_ ,, „, >'(v^i)(T,x'-(t-x)u',ll')dT 

,, // / , // Its -vlv")it-tn)-l{ i^(^/7li)(-r,x"-(t-T)u",i.")£iT| 

-\-to<t^^{x".^")h(to,x" ~(t-to)v",v")e °' J'o '■■-"' u, '^ ' | 

I if r I^ J, I T^ f ' '' M ^ ' /* > ' '\ -i-Ci-'Xt-s)-/' ,y(y7r-^)(T,a;'-(t-T)i.',i.')dT , 

+ / {KTj;h + wT^{ — , — )\{s,x —{t — s}v,v}e "• " ' Js vv-^^'^ ^ ' ' ' ds 

iymax{0,t-to-tb(=^'."')} ™ ^ 

f* r r^ 1, I T^ / '' '^ M / " ^i ^ " "\ -■'("")(*- = )-/.* ''(VJ» — )(T,x"-(t-T)i.",i.")dT J 

— / {A^o/i + ii^r^^ — , — )}{s.x ~ [t — s)v ,v }e ^ ' ^s »vf-^i,'\ V ^ ' ^ ^ 

imax{0,t-to-tb(="".i'")} ™ ™ 

It is easy to see that if t — io > ib(a;', v') then as (5 — >■ we have 

t — tToix! , v') -^ to , Xb(x', w') — ^ Xq , 

and if t — to < t\,{x' , v') then as (5 —> we have 

x' — {t — to)v' — > Xq. 

Therefore the first four Unes converge to \h]t^,xo,vo x e-''(''«)(*-*'')"-''4 ''(^™)(^^^«-(*«-^)"«^''«)'*^. For the last two lines, 
using the continuity of Kyjh,T{ — , — ), i/{y/Jl—) we conclude that it converges to zero. Therefore we have 

\h(iV <r [?,! ^-'^(i'n)(t-to)-/t' i'(\/M^)(T,a:o-(to-T)t)o,fo)(iT 

< \h] X p-(c7;7-c™c"ll''olU)(i+l''ol)^(t-to) 



where we used 



i/w{v) = / / B{v — u,uj)e 'I- w ^{u)dujdu (113) 



with ^{i + \v\y<„^{v)<cui + \v\r- (114) 



Remark that Proof of ( 18 ) is valid for in-fiow, diffuse and bounce-back cases. 



Proof of ( [20| ) 

Assume [h(to)]xo,vo 7^ and ip G (O,ib(a^0j ~^o)) with (a;o,uo) G 7o- Further assume that the boundary 951 is strictly 

concave at xo along vo direction ( [l9| . 

Step 1 Claim : We can choose sequences {t'^,x'^,v'„), {t'^,x'^,v'^) G [0, oo) x fj x R-^ n B{{to,xo,vo); ^)\{to,xo,vo) such 

that lim„^o, |/i(4,<,<) - /i«,x;;,OI > '2iHto)]xo,vo ^ 0. 

From [h{to)]x„.va 7^ we may assume 

3 

sup \h{to,Xo,Vo) - h{to,Xo,Vo)\ > -7[h{to)]xo,vo 7^ 0, (115) 

i3:'o,v'o),ix'oyQ)eB{{xQ,vo);^)\{xa,vo) 

for all n £ N. And for each n G N we can choose desired sequences. 
Step 2 Claim : For given e > 0, up to subsequence we may assume that 

{xb{x[„v'J,vl^) e B{{xo,vo);e)\& , {xb{x';„v';),v'^) (^ B{{xo,vo);e)U <3 for all n G N. (116) 

We remark that a continuity G(t, x, v) — w{v)g{t, x, v) on [0, oo) x {j_ U 7f }, i.e. 

[G\[o,ac,)x^Jto,xo.vo^w{vo)[9\[o,ryo)xyJto.xo.vo^O for aU (io,a;o,wo) G [0,oo) X {7_ U7o} (117) 



is crucially used in this step. In order to show the final goal (116) of this step, we need to prove following statement. 
Assume {xq, uq) G 7^ and t\y{xo, vq) > to. Then for sufficiently small £ > there exists iV > such that 

if (x,v) G B((xo,vq); — ) for n > N and x-^ix, v) 4. B((xo, vq); e) then we have tiy(x, v) > to- (118) 

n 
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We will prove (1181 later and show (116) using (118). It suffices to show that there are only finite n G N such that 

ixUx'n,v'J,v:^) e B{{xo,vo);-)\e , (a;b«,0,0 e i?((a;o,i;o); -)\© (119) 

n n 

or {xk{x'^,v'J,v'JiB{{xo,vo);^)ii& , (a^bK, <), O ^ S((xo,i)o); ^) U 0. (120) 



Suppose there are infinitely many n' e N satisfying (119). If e > is sufficiently small then (119) implies that 



dr 



to > tbi^'j^, , v!^,) and Iq > tb{x'^,,v",). The Boltzmann solution h at {to,x'^,,v'^,) is 
/i(to,x„,,t;„,) = n-(^o-^b(a;„,,'y„0>2;b(a;„',w„')>«n')e "„,.„,; 

and a similar representation for h{tQ^x'^,^v'^^,). Compare representations of h{ta^x'^,^v'^,) and h{tf) , x'^^, , v'^,) to conclude 

lim |/i(io,<',<')-/i(io, <',<')! == lim \Kta-ti,{x'^,,v'^,),X],{x'^,,v'^^,),v'^,)-h{tQ-txy{xl,,v'^,),Xb{x'^,,v'^,),v'^,)\ 

n'— >-C30 n'— >-C30 



< [^l[0. 



oo)X7_Jto-tb(a;o,«o),a:b(a;o,''o),i'o "^ ^ 



i'(fo)(to-«b(2:o,i'o))-/t'p''_ti_^(^jj_„p) i'(yM^)(-r,a;o-(to-T-)t)o,-uo)dT 



where we used the continuity of v{,JJi—^) and r+( — , — ). Further using the in-flow boundary condition h\^_ = wg, we 
have 

lim \Kt,,x'^,,v',,)~h{t,,xl,X[.)\<[9k,o.).,^t...o.M^^^^^^ 

where we used the continuity of g on [0, oo) x {7_ U 70}, ( |117 ) at the last line. This is contradict because we choose the 
sequences {x'^,,v'^,), {x'^,,v'^,) satisfying lim„_^oo \h{to,x'„,,v'^,) - h{to,x'^,,v'J^,)\ >h\Hto)]xo,vo ?^ in Step 1. 
Now suppose there are infinitely many n" G N satisfying (120). Because of (118) we have to < tb(a;^//, w^,/) and 
to < tb(x^„,w^'„). The Boltzmann solution h at {to , x'^^,, , v'^^,,) is 

h{to,x'„,,,v'„,,) = ho{x'„,,-tov'„,,,v'„,,),v'^,,)e~''^<'''>'o~f^°HV^^^^^^^^^ 

+ /*^if„/l + r+(^,^)}(s,x;,, -(to-s)<.,<,Oe-''K")(*o-«)-/^Kv/Al^)(r,<,,-(to-rK,,,.;„)dr^^^ 
Jo WW 

and same representation for /i(to,x"„, w^'„). Using the continuity of ho we see that 
lim \h{to,x'^„,v'^„) -h{to,x'^n,v'l„)\ 

n— >C30 

= hxa. \ho{x^„ -tov^„,v^„)-ho{x^„ -tov^„,v^„)\w{vo)e Jto-'b(-o.-o) 

n—>oo 

= 0, 
which is also contradiction. 



dr 



Now we prove (118). We can choose e > sufficiently small so that dVt n B{xo;£) = {{xi,X2j^{xi,X2)) G B{xo]e)}. 
From tb(xo, Vq) > to we know that a line segment between xq and xo — toXo has only one intersection point xo with dfl, 
i.e.xo, Xo — tovo n dV, = {xo}. Furthermore we can choose g > so large that IJsefo t 1 ^(^0 — svo] g) n dil C B{xo; e). 
Choose TV e N sufficiently large so that x,x — tov C UsgIo tol ^(-^o — svo] g) for all [x, v) e B{{xo, vq); ^). If Xb{x, v) ^ 
B{{xo, vo)',£) then x,x — tov H dfl = and this implies tb(2;, v) > to- 

Step 3 Claim : Choose t > so that t — to G [0, tb(xo, — Wo)) and denote x = Xq + (t — to)vo, v ~ Vo- Then 
there exists iV S N so that t — to < tb(a;^, —v'^) for all n> N . 

Using (118 1 we only have to prove Xb{x'^, —v'.^) ^ i3((xo, — wo);e)- From ( |116 ) we know that Xb(a:;^,w^) G B{xo]£)- We 
assume that VL(^ B{xo]£) = {a; G B{xo]£) : x^ > $(a;i,X2)} and n{xo) — (0,0,-1) and vo ~ |wo|(l,0,G). Let's define 

^(s) = mx'„)i + S«)i, {x'„)2 + S«)2) - («)3 + S«)2). 

Since x^ G fJ we have *(0) < and ^(tb(a;^, — u^)) = = ^(— tb(a;^, u^J). Because of the strict concavity along vq 
direction at xo (19), for sufficiently large n so that {Xj^,v'„) '^ (xo,wo) we have 

■ «)i ■ 



* (S) = ((«jl,(02)( qJ^^^ q2^^ 



{V'nh 



<-^Cx 
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where the Hessian of $ is evaluated at {{x',^)i + s{v'^)i,{x',^)2+s{v',^)2)- Since {x'^+sv'j^ : s € {—ti,{x'j^,v'^),tb{x'j^,—v',^))} C 
n we have *(s) < for s e {-tb{x'^,v'J,tb{x',„-v'J). Therefore ^'{-tb{x'^,v'J) < and ^'{tb{x',^,-v'J) > 0. This is 
contradiction because 

The consequence of this step is that for n > N we have a representation of h at (i, x, v) 

hit,x'^ + it-toH,,v'J = Mio,<,<)e-'^('^")(*-*°)--^o''(V7^^)(-'<+(-*'')'^"'<)''- (121) 



(122) 



Step 4 Claim : For given e > there exists 6 > so that if \{y,u) — (xo,wo)| < ^ and |(a;, u) — (a;o,fo)| < 5 and 
^0 < tbiy,u) and to < ^b(2;,f) then 

[^(icy,-")- /i(io,a;,w)| < £• 
We have h{to,y,u) = /io(y - tou,u)e-''(«)*o-/o° ''(v^S)Ky-(to-r)u,«)dr 

+ riK^h + r+( ^ , ^ )}(S, y - (io - S)U, u)e-^("H*o-.)-r° Kyj^^)(r,l;-(to-r)n,«)dr^^^ 

Jo WW 

and similarly h{to,x,v) = ho{x - tow, w)e-''(")*«^-''o° ''(v^^)Kx-(to-r)«,t;)dr 

+ riK^h + r+(-, ^ )}(S, X - (to - s)«,«)e-''('')(*"-^)-/='" .(VM^)(r,.-(to-r).,.)dr^^_ 

Jo w; w 

Let's compare the arguments of two representations : 

\{y-tou,u) - {x-tov,v)\ <2{l + tQ)S for ho, 

|(r, y- (io -r)u,w) - (r,a;- (io -t)u,u)| < 2(1 + fo)^ for v{^—), 

w 

\{s,y-{tQ-s)u,u)-{s,x~{to~s)v,v)\<2{l + to)S for K^h + r+{-,-). 

w w 



Using the continuity of h^, i^{y/Jl-^), Kyjh and r+(;^, ^) we can choose desired e > to conclude (122). 

Step 5 Claim : Choose t > so that t € [toj^o + ^bla^Oi — "^o)) and denote x = xq + {t ~ to)vQ, v = vq. Let 
£ < jQ[h{to)]xo.vo and (5 > be chosen in Step 4. Then we can choose m" e fi so that |u" — -y"] < (5 and to < fb(a;", u") 
and t ~ to < th{x'^, —u'^)- 

If there are infinitely many u'^ so that to < tb(a;^, wj() and t — to < tb(a;", — u") then up to subsequence we can define 
"n — '^n- Therefore we may assume t — to > tb(a:", —v'j[) for all n G N. We assume that f2 n i?(xo; e) — {x € B{xo;s) : 
X3 > ^{xi,X2)} and n{xo) — (0,0,-1) and vo = |t;o|(l,0,0). Now we illustrate how to choose such a u'^. Denote 
x'^ — X = (xi, X2, X3) and w" = (wi, f2, vs)- First wc will choose (ui, U2, U3) and s > so that 



n(xb(x, — u)) ■ u — 0, 



(123) 



and a:b(a;, — w) = (xi + s 



Us d 



X2 



^Jul+u: 



^,$(xi +S 



\/"?+"i 



,a;2 



Ml 



, - -<^{xi+s- , 

.Jul + ul ds ^u\ + u\ 



X2 



U2 



^(xx 



= )). The condition (123) implies that 



-\/"?+^ 



r)-a;3 



\A^F^ 



(124) 



In order to use the implicit function theorem we define 



^{xi,X2,X'i\Ux,U2;s) = $(xi 



\f<^' 



,X2 



U2 



\f<^' 



X3 



Ul 



^MV' 



= d,j:^<^{xi + S 



Ul 



\/u?+' 



, a;2 + s 



U2 



r) 



U2 



and compute , using ( 19 ) 



< 



Ul 



"2 



\/u? + ul V'wi~+~ 



^2 (i) 9 r9 $ 



f5^2*(a;i 



Ul 



y^^ 



f,X2 + S 



U2 



xA^ 



Vm?+„2' ^ul+ul' V 5-l5:.2^ 92^* 






< 



:(-'x 



(125) 
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for X ^ xq, V ^ vq and the Hessian is evaluated at (xi + s , "^ o ^,a^2 + s—!==^). Hence s = s(xi,X2,X3:wi,'W2) is a 

smooth function near a; ~ a;o and (mi,W2) ^ (^1,^2)- In order to study the behavior of s we use the Taylor's expansion 
: from ^(xi, X2, ^3; ui,U2', s) = we have 



Hxi,X2)-X3^ J (m,U2) ( A ^/^ rh ^^^2''' 



-uf + U2 



'^Xi'-'X2^ '-'X2 



Ul 
U2 






Ml 
U2 



s\ 



(*) 



(*=^) 



where the Hessian (*) is evaluated at [xi + s*^==^,X2 + ■s* / "^ 1) and the Hessian (**) is evaluated at [xi 



■ \/«f+«| 



■ -y/«f+i 



, a;2 + s* 



of the above equation converges to 



) with s^,s^,^ € (0,s). For x ^ xq and (wi,U2) ~ (^i5i'2) we know that the right hand side 



2{v 



Hence we have a control of s, i.e 



J / , „ ^ f 5^1 *((2;o)i, (2^0)2) d^,d^,<^{{xo)i,ixo)2) \ f n \ ,f. 

2+^2)i^l'"2J (^ a,,a,,$((xo)l,(xo)2) 92^$((xo)l, (0:0)2) j U2 j ^"• 



c 



\<^{xi,X2)-X3\2 <S< C\^{xi,X2) ~X3\- 



(126) 



From (1241, U3 = ^ul + ul£^{xi 



s ,t' „ ,X2 + s ,l\ J equals 



%/«?+"! ' %/«?+'' 



2 I 21 •\/mi+^ 

^1+^2 I M2 



^f4 + ~^ 



a.,<i>(xi,a:2) + ^=|i=^a2^$(x'i,x^)s+ ^=S^a,,9,,$(a;'i,4)s 

a.,<i>(xi, 0:2) + ^f=,a,,a,,<i>(x'i, 4)s + ^^£=^92 $(^, , ^^ 



(127) 



where x'^ ~ xi + 3^ , "^ =, x'2 — X2 + s' j "''^ = for some < s' < s < C|$(a;i, 2:2) — xsj 2 . Using the smoothness of <!> 

V'"l''""2 



we can bound (1271 as 



y^nf+wf 



^|(ui,U2)| (|(a;i,a;2)| + |$(xi,a;2) - xal') < (1271 < C\(u^,U2)\ (|(xi, 0:2)1 + |$(a;i,a;2) - a^gl^) . (128) 

To sum for fixed x and direction . "'" (u\^U2) we can choose U3 such that ri(xb(x, — (ui, ^2,1*3))) • (ui,U2,M3) = 



111 -\-ii 

and U3 is controlled by (128). Finally we choose (1*1,1*2) — \ \, I (fi, ^^2) and find the corresponding M3 so that 
|w| — \(u\,U2,U'i)\. Define mJ( = — w + 2(i; • |("^'"'''"^)| )('»i;'»2, "a)- Then we have desired u" for sufficiently large n G N. 

Step 6 To sum for (i,x^ + (^ — ^o)Mnj'"n) we have t — to < tiy{x'^, —u'^) and io < tbixn,u'^) and |/i(to, a^^jwJO — 
h{to,x'^,v'J^)\ < YQ[h{tQ)]^g yg. Hence the representation of the Boltzmann solution h at {t,x'^ + {t — to)u'n,u'^) is given 

by 

hit, x^ + it- toK, - h{t„, <, Oe-''^"")^*^*"'^^^ -(vj^^)(-'<+(-^ow:.u':)dr 



h h , 

w' w' 



+ / {i^„/l + wr{-, -)}{S, x:: + (S - toX, Oe^''^)^*-^)-/^- KyjI^)(r.<' + (r-to)<,<)dr^5. 



Using (121) we have 



)dT 



lim |/i(t, x'^ + {t- to)v'^,v'J ~ h{t, x';, + {t- ioX, 01 

n— J-oo 

= lim \h{to,x'„,v'J - /j(tp^^'^,u^^)|e-''("«)(*-*'')-/*>(VMS)K.o+(r-to).o,.o) 

n— J-00 

> { lim \h{to,x'^,v'J-hito,x:,v'::)\- hm |Mio,<,<)-Mio,x::,<)|}e-''("'')(*-*'')--^>''(^^)(^'^°+(^ 

— A ["■\''Q)\xo,vof^ " J 



-ta)va,VQ)dT 



which implies that 



[h(t)]x,. > \[h{to)W.„ X e-(c.+c'c„||.olU)(i+|.|)-(*-*o) ^ 0. 



Remark Through Step 1 to Step 6, we only used the in-flow boundary datum g explicitely in Step 2. All the other steps 
are valid for diffuse and bounce-back boundary condition cases. In Step 2, we only used (111) the continuity of G = w g 
on [0,00) X {7_ U 7^ }. Therefore, if we can show the continuity of F on [0,cxi) x {7_ U 7q } then we can prove (20). 
For diffuse and bounce-back boundary we will prove such a continuity to conclude (20). 
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In this section, we consider the Unearized Boltzniann equation (911 with the diffuse boundary condition (93). In spite 



of the averaging effect of diffuse reflection operator, we can observe the formation and propagation of discontinuity. 
Continuity away from 2) is also established. 



5.1 Formation of Discontinuity 

We prove Part 2 of Theorcmll] The idea of proof is similar to in-flow case but we also use |wo| not only to as a parameter. 
Without loss of generality we may assume xq — (0, 0, 0) and vq ~ (|wo|, 0, 0) and {xq, vo) S 7o • Locally the boundary is 
a graph, i.e. nnB{0;6) = {(xi, Xa, 2:3) G B{0;S) : X3 > <^{xi,X2)} and $(^,0) < for f e (-(5,J)\{0}. (See Figure 3) 



Assume that ||/io||oo < ^ is sufficiently small so that the global solution h of (91 1 with diffuse boundary (93) has 



a uniform bound ( [95| , from Theorem 4 of [IH]. Choose to G (0,min{(5, tbla^o, ~t'o)}) sufficiently small and \vq\ > 
sufficiently large so that 



-< f e-''(l"«l)*°-toCkC"-(l-e-''(l'"'l)*'')Cr(C")^-C"— ^, , 

2 \ w(vo) J{v'^>Q} 



w{v')da{v') , 



(129) 



where i^{\v\) = i/{v) and Ck and Cr from (33) and (34|. More precisely, first choose \vo\ > large enough to have 

1 



then choose to > as 

< to = min 



w{vo) 
S tb(xo,-wo) S 



{i+p'\vo\r ^ 1 



e * 



IOC" 



1 , / 10Cr(C")^ 
log 



|i;o| ' lO^d^^ol) ' lOCkC" i^d^ol) " yiOCriC')^ - 1 



Assume the condition for initial datum ho : there is sufficiently small 6' = S'{il, to|wo|) > such that i?((— to|wo|, 0, 0), S') C 
fl and 

/io(xo,i'o) = ||/io||oo>0 for {x,v) e Bii~tQ\volO,0);S') X B{{\vo\,0,0);S'). (130) 

We claim that the Boltzmann solution h with such initial datum Hq is not continuous at (to, xo, vq) = (to, (0, 0, 0), {\vo\, 0, 0)). 
We will use a contradiction argument : Assume the Boltzmann solution h is continuous at (to,a;o,wo), i.e ( |98[ ) is valid. 
Choose sequences of points {x'„,v'J = ((0,0, ^), (|wo|, 0, 0)) and (a;„,u„) = ((i, 0, $(^,0)), -^i^(|i;o|,0, -^0)^ Because 

of our choice, for sufficiently large n G N, we have 

(x;-to<,<) = ((-to|i;o|,0,-),(|t^o|,0,0)) e B{{-to\volO,0);d') x B((|z;o|,0,0);<5')- 



Hence the Boltzmann solution at (to,a;^,f^) is 



h h. 



h{to,x'„,v'„) = ho{x'^~tov:^,v'Je-''^<^'°+ e-'^<^^''>-^U K^h + wT{-, -)\ iT,x'„ ~ v{to ~ T),v'„)d7 

Jo I WW) 

= \\ho\\ooe--^^<^^'" + [ ° e-''(l<J)(*«--) U^h + wTi-, -)] (r, x - <(to - r),v'JdT. 
./o w w \ 



Using the diffuse boundary condition (93 1, the Boltzmann solution at (to,x„,w„) € [0, 00) x 7_ is 

1 



h{to,Xn,Vn) 



*(|fo|) Jv{x„) 



h{to, Xn, v')w{v')da(y'). 



Using a pointwise boundedness (95 1 of h, and ||ft.o||oo < 1, we can estimate 

\h{to,x'.^,v'n) -/i(to,a;„,w„)| 

> I ll/iolUe-'^d''"!)*''- /''{CkC'||/io||oo + ^«)e-''(''")(*»-"^Cr(C")'||/io||L}rf^-C^'ll'io||oo^7^ f w{v')da{v')\ 
Jo w{\vo\) Jv 



> ||/io||ooe-''(l"«l'*« - toCkC'll/iolloo - (1 - e"'^(l''«l)*«)Cr(C')'||/io||L " C'\\ho\ 



1 



w'lPolJ Jv 



w{v')da{v') 



ll/iolloo fe-'^d'^"!)*" - toCkC - (1 - e-''(l'^''l)*°)Cr(C')' - C— ^ / w{v')da{v')] > 
\ w{\vo\) Jv J 



l^ollc 



^0, 



which is contradiction to (98). 
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5.2 Continuity away from 2) 

Instead of using the argument of |13) to show continuity in the case of diffuse reflection boundary condition we will use 

the sequence (103) with the boundary condition (131) and Lemma 11 This argument also gives a new proof of the 

continuity of Boltzmann solution in a strictly convex domain with simpler way than [13j . 

Proof of 2 of Theorem [s] 

We will use the sequence (1031 with h"^^^\t=o = ho with following boundary condition 

1 



h"'+'{t,x,v) 



Wiv) Jv{x) 



h"'{t,x,v')w{v')dcT, 



with (i,x, w) S 7_. 
Step 1 : We claim that 



w{v) 



V{x) 



h"^{t, X, v')w{v')da{v'), 



(131) 



(132) 



is continuous function on [0,T] x 7 even if ft," G L°°{[0,T] x fix R^) is only continuous on [0,T] x fix M^\0. We will 
show as {t, x,v) — >■ {t,x, v), 



1 



wiv) Jvi^^ 



K^{t,x,v')w{v')da{v') 



w{v) 



V(2) 



h'^{t,x,v')w{v')da{v'). 



(133) 



Using the fact |V(a;)\V(x)|, |V(x)\V(x)| — > as x — >■ a; and the exponentially decay weight function of wda it suffices to 
show that 

/ {w{vy^h'"'{t,x,v')w{v')-w{v)-^h"'{t,x,v')w{v')}da{v'), (134) 

JV(x)nV(x)n{\v'\<M} 

for sufficiently large M > 0. Using Lemma^we can choose open set Ux C {v' e R? : \v'\ < M} so that \Ux\ is small and 

/i™ is uniformly continuous on {|i;'| < M}\Ux- Therefore we can make S\;(x')r\v(x')nUv'\<M\nu '^niall using the smallness 

of Ux and make /v(-2;)nv(2')nllD'l<M>\(7 small using the uniformly continuity of ft™ on {|t;'| < M}\Ux- Hence (1331 is 

valid. 

Step 2 : We claim 

ft.* is a continuous function in C^ (135) 



for alH G N where €.t is defined in (106 1. By induction choose hP — Q and (135) is satisfied for i = 0. Assume (135) for 
ah i = 0,1,2,...,™. Let wr_ f^,^^^^^ =:i/(^)ft™+i. Then the equation of ft™+i is 

/ft'"\ /ft™ ft'" 

{dt + vyx + v{v) + i^[—] }ft™+' = if^ft™ + w;r+ — , — 
\ w J \ w w 

From Theorem W and Corollary p| we know that v (^) and wT^ (^, ^) are both continuous in [0,T] x il x M.^. 
Because of Step 1 we know that :i^jy) Jyr^) h'^{t,x,v')w{v')da{v') is also continuous function on [0,r] x 7. So we can 
apply Lemma [TT] to conclude ( |135[ ) is valid for i = ?ti + 1. 
Step 3 : We claim ft'" is a Cauchy sequence in €-t for some small T > 0. First we will compute some constants 



explicitly. From ( 6 ) the normalized constant c^ is 
we can compute tne right hand side of above term : 



) 2 /"OO 2 /"OO 

dvi Vie~^ I dv2 e~^ I dv^^ e~^ 
J —00 



Therefore we have c,, — 7^. Next we will show 



l7i(x)-v'>o^ 2 {n{x) ■ v'}dv' . Choose n(x) = (1,0,0) and then 

{ -e'^ ] dvi X {V2tt)'^ ^2tt -e' ^ 

dvi \ J 



2tt. 



W{v) Jy'.n(x)>0 



w{v')da{v') < Cpp^^~\ 



(136) 



where w{v) ^ = (1 + p'^\v\'^)^e ^ . We follow the computation of Lemma 25 in [13j. For ^jt^, in the case of /3p^ > j 
we can see that ■w{v)~^ has a maximum value at |w| 



- . / M^ which is 



il + p'\v\re 



'\v\-- 



ipp^-l 



'e-^e^p^^, 



(137) 
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* < 



'w{v')da{v') — —- 

v'-n{x)>0 ^^ 

1 

2^ 



.'12 |„'|: 



(l + p2|i,'|2)-^e^^e-^^«W^^' 



'u'i>0 



Cpp 



ui>0 
4 



27r7„,>o(l + |ti|2)^-i 



du 



where /? > 2 and combining with (|137|) we conclude (|136|). 
First we will show a boundf 

Lemma 12 Lei ft,™ 6e a solutu 
SMc/i i/ia< «/ ||/io||oo < S then 



First we will show a boundedness ( |108 ). 
Lemma 12 Lei ft.™ 6e a solution of (103) with ft.7=o — ''-o o,i^d the boundary condition (131). Then there exist T^,,C,S > 



sup ||ft™(s)||oo < C||/io||oo for all to e N. 

0<s<T, 

Proof. We will use mathematical induction. Choose ft° — Hq and assume ||/io||oo < ^ and 

sup \\h\s)\\oo < C'WhoWoo, (138) 

0<s<T. 

for i — 0, 1, 2, ...,TO,, where 6,C,T^ > will be determined later. From Lemma 24 of [13 the representation of /i™+^ 
which is a solution of ( 103 ) with the boundary condition ( 131 1 is given by 

ft 



h"'+\t,x,v) ^ lt,<oit,x,v)i.ho{x~tv,v)e"''^''^'+ e"''^''^^'-'U"\s,x ^ {t - s)v,v)ds\ 



(139) 



[initial data] 



-i/(t))(t-ti) 



+ lo<tAt,^,v){ e-''^-'>^'-''>q"'{s,x-{t-s)v,v)ds+- 



H 



}, (140) 



II 



where g™ was defined ( 104 ) and 

fc 



H = y2Ui+i<o<tiho{xi -tivi,vi)dT,i{0) + V] / lt,^i<o<t,g'""'(s, a;; - {ti - s)vi,vi)dY,i{s)ds 
1=1 ' ^ ' 1=1 io . 



[initial data] 



III 



fc /-t, 



y^ / lo<t,+ig™ '^{s,xi- {ti- s)vi,vi)dT.i{s)ds + lo<tk+ih"^ '''^^{tk+i,Xk+i,Vk)dT.k{tk+i) ■ 



(141) 



(142) 



j^i yii+i 



IV 



[many bounces] 



Here d^kitk+i) is evaluated at s = t^+i of 

dS,(5) = {nf^;+idaj}{e-''('")(*'-^)iI;(«i)dc7i}n;.-\{e-'^(^^)(*-*^+^)da,}. 



First we can estimate [initial data] in (139) and (141), 



/ < lt,<o\ho{x - tv, v)\ + — — Y^ lt,^^<o<t, |fto( 

< maxi 1, -7-- max / w(w;)iicri...(i(7fc ^||/io||oc 

I w(w) i</<fe7n,"=iV, J 



X; - tivi,vi)\w{vi) > dui.-.dak 



< 

where we used ( 136 1 . 

Next we estimate [many bounces] term in ( 142 ) which is crucial estimate in this proof. We use Lemma 23 in |13| to 

bound a contribution of [many bounces] term in ( 142 1 in the last term of ( 140 ) by 

-^^T^ / ^{t^+i{t,x,v,vuV2,...,vk)>o}wivk)dakdak-i--.d(7i x sup ||ft™"''+^(s)||oo 

w[v) Jyi'^^^ Vj Q<s<t 

- ~^r~^ / w{vk)dak / 'i-{tkit,x,v,vi,...,vk-i)>o}d(Tk-i---dai x sup ||/i'""''^^(s)||oo 

'w[v) Jy J]!''!} V, 0<s<t 



'mzt 



< Cpp 



2/3-4 



C^P 



5/4 



sup W-''+\s)\\^<Cpp- 

Q<s<t 



2^-4 



C2P" 



C\\ho\\c 
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where we used (136) at the last step. The remainders I, II, III and IV arc contributions of q™ , . . . , q™ '^. We introduce 
a notation 



n, = tCk sup \\h\s)\\^+Cr sup W{s)\\^( sup \\h\s)\\^ + sup W+\s) 

a<s<t 0<s<t \0<s<t 0<s<t 

< C||;io|U(CkT, +2CCr||/iolU), 
where the above inequahty holds for < t < T* and i — 0,1,2, ...,m — 1 and 

n^<{T.Ci, + CrC\\ho\\oo)C\\ho\\oo + CrC\\ho\U sup W+\s)\U 

0<s<T, 



(143) 
(144) 

(145) 



where we used the induction hypothesis (138) for (144) and (145). Easily we have 



I, II < Hm, 



III 



IV < V] -r-r^ / dai.. / dai-i / dai+i.. / dat / / H,, 

l^^ W(WJ Jy^ Jy^_^ Jy^^^ Jy^ Jy^ Jq 

< y]^™-'^ / wivi)dcTi<Cpp'"'-^y2^rn-l- 



-le 



-v(vi){ti~s) ~ 



'w{vi)dsdai 



To summarize, we can estimate all terms of representation of h™-^^{t, x, v) in (139) to obtain 



\h"'+\t,x,v)\ < \\ho\\oo{c[2T,Ck + 2CrC\\ho\\oo + kCfip^^-^iCi,T,+2CCr\\ho\\oo) + Cpp''^-^ 

+l + Cpp''^'-^} + CrC\\ho\\oo sup ||/i"+i(s)|U. 







Csp-^ 



0<s<T, 



Choose k = p^/*. Choose p > sufficiently large so that Cpp^^ ^ {W ^ — ^ ^'^'^ then choose T* > sufficiently 
small so that T* x Cr(l + C,ip^/'^p^^-'^) < ^ and then choose C > sufficiently large C > 10(1 + (7/3^^/3-4) a^j choose 
(5 = min j ^^^ ^ , gg|=^ (CC^p^/^/?^^^^) >. Finally assume ||ft.o||oo < <5- Then we have 



1 1 ; 71 

sup 1 1 h 

0<s<T, 



"'(5)||co < 1 ^ /.n. II l|feol|oo|l+' 

i - c-rO||ftol|oo ^ 



Cpp^P-^ + C 



Cpp'^-'i-A^'^" +tCr(l + C,p^/V^-^) 



20 
- 19"''"- 110 



00 ^ tt: ~^ C 



+CrC||/io||oo + 2CrCpp''/^p^^-*C\\ho\\oc] } 

<C\\ho\\oc- 



1111 
30 ^30 ^20 ^15 



Next we will show that /i™ is a Cauchy sequence in L° 



Lemma 13 Let /i™ be a solution of (103) with h™'^^\t^Q = /iq and the boundary condition (131). Then there exist 
T^,C,S> so that j/||/io||oo < S then /i™ is Cauchy m L°°{[0,T^] xClx R^). 

Proof. The equation of /i'"+^ — /i™ is 

with {/i"+i-/i"}L=o = , {h'^+^-h'^jL^^-^ [ {h'^(t,x,v')-h"'-\t,x,v')}w(v')dc7(v'), 

Wiv) J(^) 



where g™ is defined at (110 1. From Lemma 24 of [13] we have the representation 

{h"'+^ ^h'^}{t,x,v) ^ lt,<oit,x,v) [ e-''^''^^'-'''^q"'{s,x-it-s)v,v)ds 

Jo 



(146) 



+^o<tiit,x,v)\ / e 



^uiv){t-s)~m(^g^ a; - (t - s)v, v)ds +' 



-iy{v)(t-ti) 



H 



nu V, 



}■ 



II 
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where 

^ = V] / 'yti+i<o<tiq"^^\s,xi~ {ti- s)vi,vi)(n:i{s)ds 

1=1 il ' 
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III 



fe /.*, 



Y, / lo<t,^,qr'\s,xi - {ti - s)vi,vi)d^i{s)ds + U<t,^^{K^-''+^ - h^-''}{tk+i,Xk+i,Vk)dEkitk+i) . 



/=iy+i 



[[many bounces]] 



IV 



First using Lemma 24 of [T3|, we estimate [[many bounces]] term for sufficiently large A; > by 
f 



'w{v) 

< Gpp'^'U-} sup \\{h 

L^J o<s<t 



{t^.(t,X,V,Vi,...,Vk 



_^)>a}dak-i-dai x sup \\{h 



rn — /c+l im—k 



~f . I -{t,+i(t.x.v,vuV2,...,vk)>o}w{vk)dakdak-i.-d(TiX sup ||{ft.™ + - /i" }(s 

^(u) Jn^^i Vj ' o<s<t 

— I w{vk)dak I 



0<s<t 



/i™-n(s)l|c 



Easily we have I, fl < 5^™, III, IV < Cap^'^-'^dnm-i where 



6n, = tC^ sup ||{/i'-/i^-i}(s)||oo+C||/io||ooCr( sup \\{h' - h'-^}{s)\\^ + sup ||{/i'+' - /^*}(s)|U) 

0<s<t 0<s<t 0<s<t 



< - sup ||{/i*-/i-^}(s)|U+ sup ||{/z»+^-/l*}(s)||, 
4 Lo<s<t Q<s<t 

with r = 4max{iCk,C||/io||ooC'r}. 

To summarize, we can estimate all terms of representation of h"^^^(t, x, v) — h™{t, x, v) in ( 146 1 for any m > k to obtain 



sup \\{h^+^ - h"^}is)\\^ < -^i^Cpp'^-^J2i ™P ll{^ 

0<s<t l — ZTlZ J^^O<s<t 



m—L urn- 



-l-l 



}( 



s)lloo + sup \\{h 

Q<s<t 



m— Z+1 im—l 



h-'-'mWo.) 



- sup ||{/i™-/i™-l}(,s) 
^ 0<s<t 



C, 



PP 



2^-4 



sup ||{/i"-'=+i-/i'"-'=}( 



. 2 J o<5<t 

which is our starting point. Fix a small number r > chosen later. Choose p > sufficiently large so that 



2Cpp^'^-^ {1} ^^ < I and then choose t > so small that ^C/jp2/9-4 <- f ^^^^j _r^ ^ f rpj^^^^ ^^ j^^^^ 
sup ||{/i'"+^-/i"}(s)|U<t( sup ||{/i™-/i"-i}(s)|U + ...+ sup ||{/i"-'=+i-/i"-fc}(s)||, 

0<s<t l0<s<t 0<s<t 



(147) 



Using (147) for m, j G N so that m — {i + l)fc > and j = 0, 1, ..., m — 1 it is easy to show 



sup ||{/i™-'fe+i+J - h"'-'''+^}{s)\\oo < 

0<s<t 



(l + fyl sup ||{/i"-*'=-/i"-''=-i}(s)|U + ...+ sup ||{/i™-(»+i)'=+i-/i™-('+i)fc}(, 



0<s<t 



We apply the above inequality term by term in ( 147[) to have 



sup ||{/i™+l-;i™}(s)|U<T{(l+f)^--l}{sup ||{/l'"-'=-/l"-'=-l}(s)|U + ...+ sup ||{/l'"-^'=+l-/l"-^n(s)||oo} 
0<s<t 0<s<t 0<5<t 

<f{(l+f)*-l}H sup ||{/l"-'''^-/l"-'''^-l}(s)||oo + ...+ sup ||{/l™-(*+l)'= + l-/l™-('+l)'=}(s)||oo}. 
0<s<t 0<s<t 



Now we estimate 

sup \\{h"' - h"}{s)\\^ < J2 sup ||{/i™-'-/i™-'-i}0 



m—n—l 



0<s<t 
m—7i—l 



'm—ik — l—1 im—ik—l — 2\ 



< J2 r{(l + ff-ir{sup \\h 

m—n—l 

< J2 ^{(l+^)'-l}['''^^]"'{ sup ||/l2fc-/l2fe-l 



-/l™-"=-'-^||oo + ...+ sup ||/l'"^(^+l)'=-'-/l'"-('+l)'="'-l|U} 



0<s<t 



;=o 



0<s<t 



sup ||/ll-/l°||oo} 
0<s<t 
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m—n—l 

< f{(l + f)'=-l}[T^]-i y {{l + T)''-l}i'^"i^]-i'^]{sup \\h''''-h 

0<s<t 
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2fe u2k-l\ 



1=0 



sup \\h^ - h^Woo} 

0<s<t 



< T{il + f)" - l}lT^ 



^:^{ sup lift. 



2k u2k-l\ 



2-(i+f)^v<:<t 



sup ||ft^-/l°||oo}, 



0<s<t 



where we choose i = [ ™ ^ ^ ] — 1 so that m — {i + l)fc — I — 1 d [0,k). If f > is chosen sufhciently small so that 
(1 + f)'^ — 1 < 5 then {(1 + f)'^ - l}i^J^"^ — > as n — > oo which implies that 



sup ||{/i™-/i"}(s)|U^O, 

0<s<t 



(148) 



as TO, n — >■ oo. Thus /i™ is Cauchy in L°°. m 

Step 4 : We claim that h is continuous in €. Notice that T only depends on ||ft.o||oo and supq<j,<2i ||wg(s)||oo (Theorem 
1 of [13j). Using a unform bound of supq<j,^q2 ||ft(s)||ooj we can obtain the continuity of h for all time by repeating 
[0,r], [T, 2T], .... If the boundary 9rj does not include a line segment (|6| then every step is valid with [0,oo) x {JixM^}\S 
instead of € and [0,T] x {Cl x R^}\D instead of €t- 

5.3 Propagation of Discontinuity 
Proof of 2 of Theorem [2] 



Proof of (18 1 : The proof is exactly same as in-flow case in Section 4.3. 



Proof of (20) The proof is exactly same as the proof of in-flow case in Section 4.3 except Step 2. As we mentioned 
in Remark of Step 2, we need to show a continuity of a boundary datum on 7_ U 7f . In diffuse reflection boundary 
condition case, we need 

= [ ft'|[0,oo)X7_]t,i/,i) 



lini sup 

*^o t\t" e B{t;S} 

(j/'i v'), iv", v") e 7- n B({y, v); S)\{y, v) 



\h{t',y',v')-hit",y",v")\ 



lim sup 

^^° t',t" e B{t;S) 

{y\v'), (y" ,v") G 7_ n B((y, v); S)\{y, v) 



-^ [ hit', y', >o)w{^)da{^) - -^ / h{t", y", t,)^(o)da(v,) 

wW) Jv(y') W{v") Jv(y") 



for (y,w) G 7- U 7^ . This is already proven in section 5.2 Continuity away from T). 



6 Bounce-Back Boundary Condition 



In this section, we consider the linear Boltzmann equation (91 1 with the bounce-back boundary condition (94) 



6.1 Formation of Discontinuity 

We prove part 3 of Theorem fll Without loss of generality we may assume xq = (0,0,0) and vq = (1,0,0) and 
(a:o,i'o) S 7o. Locally the boundary is a graph, i.e. ft Ci B{0;S) = {(xi, 0:2, 2:3) G B{0-,S) : X3 > ^{xi,X2)}- The 
condition {xo,vn) G 7o implies tb(a;o,wo) 7^ and tb(a;o,— ^o) 7^ which means $(^,0) < for ^ e (— (5, (5)\{0}. (See 
Figure 3) 



Assume that ||ftol|oo < ^ is sufficiently small so that the global solution h of (91) with bounce-back boundary (94) 



has a uniform bound ( 95 ) , from Theorem 2 of |13j 
Rec 

so that 



RecaU the constants Ck and Cr from (|33|) and (|34|. Choose to e (0,min{|, bA^iz:!!ol^ tb(^o,i-o) |) sufficiently small 



1 < (^e-^'fi)*" - toCkC - (1 - e-''(i)*°)Cr(C')') 



(149) 



Assume a condition for the initial datum Iiq : there is sufficiently small 6' = 6'{fl,to) > such that 
S((-to,0,0),5'), S((to,0,0),^') C n and 

/io(a;,w) = ||fto||oo>0 for {x,v) e B((-to,0,0);6') x B{{l,0,Q);6'), 
hoix,v) = -\\ho\\oo>0 for (x,w)eB((to,0,0);<5')xB((-l,0,0);y). 



We will use a contradiction argument : Assume the Boltzmann solution h is continuous at {to,xo,vo), i.e. (98) is valid. 
Choose sequences of points {x',^,v'^) = ((0, 0, ^), (1,0,0)) and (x„,w„) = ((i, 0, $(^, 0)), ^i^(l, 0, ^)). Because of 
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our choice, for sufficiently large n S N, we have 

K - io<, <) = (Mo, 0, -), (1, 0, 0)) e B{i-to, 0, 0);S') x S((l, 0, 0); 5'), 

n 



rX 



e B((to,0,0);,5')xS((-l,0,0);(5'). 



^(-1,0,—)) 
7,2 n 



Hence the Boltzmann solution at (io,a;^,i'^) and (to,2;„,w„) is 



/i(io,<,<) = ||/iol|ooe-''«)*" + / " e-^^-^^^'^-^^K^h + wT f ^, ^')}(t,x:, - (-<)(to - r), -<)dT, 

Jo V^ "^Z 



= -WhoW^e-"^-^-')'" + / e-^(-''")(*o-^){i^„/j + «;r -, - ]}{T,Xn ~ (-i^„)(to - r),-t;„)rf 



/i /i 



Using a pointwise boundedness (951 of ft. with (33 1 and (34), we have 






2 

oo ' 
2 
oo 



Therefore using ( 149 ) 



h{to, <, v'J - h{to,Xr,,Vn) > 2||fto||oo (e-^^i)*" - ioCkC - (1 - e-'^(i)*'')Cr(C')') > ||fto||oo y^ 0, 



which is contradiction to ( 98 1 



6.2 Continuity away from l^bb 

We recall some basic facts to study the bounce-back boundary condition from '13]. 

Definition 7 JlSf (Bounce-Back Cycles) Let (t,x,v) ^ 70 U7_. Let {tQ,xo,VQ) = {t,x,v) and inductively define for 
k>l: 

{tk + l,Xk+l,Vk+l) = (tk ~ tb{xk,Vk),Xb{xk,Vk),-Vk)- 

We define the back-time cycles as: 



Xci(s; t, x,v)=Y^ l[t,^^_t^) (s){xfc + (s - tk)vk}, Vci(s; t, x, w) = ^ l[t,^^^t,) («; 



Wfe- 



Clearly, we have ffe+i = {—l)''^^v, for fc > 1, 



a;fc = 



1 - (-1)'= l + (-l)'= 



-Xl 



-X2, 



(150) 



(151) 



where xi = x — t\y{x,v)v and 3:2 = x— [2tiy{x, v) +tb(a;, — w)](— f) and let d ~ ti —t2, then t^ — ife+i = d > th{t,x,v) > 
for fc > 1, and 

ti{t,x,v) = t -ti,{x,v) , 

i2(t, X, v) = ti - tb(a;i, -yi) = ii - (tb(a;, i^) + tb(a;i, Wi)) = ti - (2tb(a;, «) + ib(a;, -v)) , 



tk+i{t,x,v) = ti - fc(2tb(a;,t;) + tb(a;, -«))■ 



(152) 



Lemma 14 /jjy Lei /iq G L°°(J7 x R'^) and <j){t,x,v) with suprg -pi^n |0(-, •, w)| < 00. There exists a unique solution 
G{t)ho of 

{dt + v-\7^+ 0}{G(i)/io} - 0, {G(0)fto} - ho, 
with the bounce-back reflection {G{t)hQ}{t, x, v) — {G(t)hQ}{t, x, —v) for x G dO.. For almost any {x, w) G Jl x R'^ \ 70, 



{G(t)fto}(i,x,«) = ^l[,,^,,,,)(0)/.o(Xei(0),K,(0))e-/o^(-^^c.(.).K,(r)M.^ 
fc 
where Xc\{t) = Xci{T;t,x,v) and Vci{t) — Vci{T;t,x,v) in (J50[). 



(153) 
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Lemma 15 (Continuity away from D^f, : Transport Equation) Let il be an open subset of M.^ with a smooth 
boundary dil and an initial datum hQ{x, v) be continuous in ilx M'^U{7_ U7+ U7q}. Also assume q{t, x, v) and (j){t, x, v) 
be continuous in the interior of [0,T] x il x M.^ and supjQ jni^nxR^ l9(^i 2;,'i;)| < cx) and swpiQrp^^^ \(/){-,-,v)\ < 00 for all 



V e 



Let h{t,x,v) be the solution of 



{dt + V ■ V:e + (j)}h — q , h{0, X, v) — ho J h\j_ It, x, v) — h{t, x, —v). 



Assume the compatibility condition on 7_ U 7q 



I- 



ho{x,v) = ho{x,-v). 



Then the Boltzmann solution h{t, x, v) is continuous on Cbb- Further, if the boundary dVl does not include a line segment 
im then h{t,x,v) is continuous on a complementary set of the discontinuity set, i.e. [0,T] x {Cl x M.'^}\'Sbb- 

Proof. The proof is similar to the proof of Lemma 16 of T?. Take any point (i, x, v) € [0, T] x f2 x M.^ and recall its 



back-time cycle and (153). Assume t^+i < < t„j. Using (1531, h(t,x,v) takes the form 






m-1 ^t 



+ E / 9(^'^fc - (^fc - ■^)Vk,Vk)e-^^-° ^*'Vi 0(r,..-(t.-r)..,.,)d-/> ^ir,.uHt.~rW,.,)dr 



k=0 •'^fe+i 



+ / q{s,Xrr,-{tm-s)v.n^,Vm)e ^-=0 Jt.+i "^^ ' ' l- ^--^ Ja 



r)!;^,^ ,v^-n)dr 



(154) 



Take any point (i, x, v) G Cbfc. By the definition of ^bb we assume that {x, v) G 57 x M'^ or (x, v) e 7_ U7q^ and we can sep- 
arate three cases : t—tb{x, v) < , {xi,{x, v), v) £ 7_U7q~ with t < 2ib(a;, v)+tt,{x, —v), and {xh{x, —v), —v) G 7_U7q^ 
with {xh{x,v),v) e 7_ U7o". 

Case of t < tb(a;, v) Simply we have h{t, x, v) = ha{x~tv, v)e- lo <P{r.x-it-T)v,v)drj^jt ^(^^ x-{t~s)v, v)eJ's Hr,x-(t-r)v,v)dT^g 
and use the continuity of q{t, x, v) and (f>{t, x, v) to conclude the continuity of h{t, x, v). 

Case of {xh{x, w), f ) G 7_ U 7q^ with t < 2th{x, v) + t^lx, —v) A representation of h{t, x, v) takes the form 



ho{xi - tivi,vi)e- ^^\^'^^''''^''^^'''^^''^~^°' '^^^''''-^''-^^''"^'^''^ + / q{s,x ~ {t - s)v,v)e- f>^^'--'^'-^>'^^''^ds 

Jti 

+ r (?(S, XI - (ti - s)vi,Vi)e- A 'i>ir,.cHt-r)v,v)dr-f^^ ^[r,x,-{t,~r)vuv,)dr ^^_ 



Thanks to Lemma 111 and Lemma l2l the condition (a;b(a^, w), w) G 7_ U 7q implies continuity of xi{x,v) 
Xb(a;, v) , ti{t, x,v) — t ~ t\:,(x, v). Therefore we can show the continuity of hit, x, v). 



Case of (a;b(a;,— t;), — w) G 7- U 7q with (a;b(a;,i;),i;) G 7- U 7q We have (154 1 for h{t,x,v). Thanks to (1511 and 

(152 1 and Lemma 111 and Lemma ^ the conditions {xb{x,—v),—v) G 7- U 7q~ and {x]3{x,v),v) G 7- U 7q~ imply 
continuity of X]^{x, v), Vk{x, v), tk{t, x, v). Therefore we can show the continuity of h{t, x, v). ■ 

Proof of Part 1 of Theorem 3 

Following the in-flow and diffuse cases, we use the iteration scheme (1031 which is equivalent to (107) with bounce-back 
boundary condition /i'"+^|^_ (i, x, t;) = h'^'^^{t,x,—v) and an initial condition /i'""'"^|t=o = ^o- 

Step 1 : We claim that h^ is a continuous function in ^bb.T for all i G N and for any T > where ^bb,T = 
£bbn{[0,T] xOxM^}. Choose /i° = and use mathematical induction. Assume h^ is continuous ^bb.T for i = 0, 1, 2, ...,7tt,. 
Apply Lemma 15 to concluse that ft,™+^ is continuous in €bb,T- 



Step 2 : We claim that there exist C > and S > such that if C||/io||oo < ^ then there exists T ~ T(C, (5) > so that 
suPo<.<r ll^"'(s)||oo < C\\ho\\oo and {/i™}^^o is Cauchy in L~([0,T] x^x E^). 

Fisrt we will show the boundedness using mathematical induction. Assume supo<5<7- ||/i'"(s)||oo < C'lJ/ioHoo where 
T > will be chosen later. Applying Lemma 14 and q correspond with v and the right hand side of ( 103 ) respectively 
to have a representation of h"^~^^{t, x, v) 



/io(^ci(0),yci(0))e-'^("'* 



-iy{v){t^s) 



{K^h'' 



wTa 



wT^ 



•m um+l 



}(s,Xci(s),Vci(s))ds, 
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where [Xci{s),Vci{s)] = [Xci{s;t,x,v),Vci{s;t,x,v)] is in (150). The above term is bounded by 



\ho\\^+tCi, sup ||/i"(s)|U + Cr sup \\h"\s)\\^ sup (||/i"(s)|U + ||/i"+'(s)||oo), 

0<s<t 0<s<t 0<s<t 



where the constants are coming from basic estimates, (33) and (34). Choose C > 4 and S < i^ and T — ^„ ^ . Then 
we have sup(,<s<T ll^™+^(s)||oo < C||/io||c 



Next we will show {/i"}^=o is Cauchy in L°° ( [0, T] x rj x M^) . Recall 9" (t, x, v) from ( [llO| ) . The equation of /i™+i - /i™ is 
(1091 with a zero initial condition (ft,™+^ "/i™)|t=o = and a bounce-back boundary condition {h"^^^ — h™-)\^ _ {t, x, v) = 



(/i™+^ - /i'")(t, a;,-w). Applying Lemma 14 to (1091 we have 



(/!'"+' -/i")(i,x,t;)- / e-''('')(*-^)g'"(s,Xci(s),t/c,(s))ds, 



where [Xci(s), Vci(s)] = [Xci(s; i, x, v), V^i(s; i, a;, v)] is in ( 150 1. Then we have exactly same estimates of in-flow case to 

conclude {ft.™} is Cauchy. 

Step 3 : Same argument as in-flow case but substitute Cbb.T , ^bb , '^bb.T , 3Dfc6 for £t , C , T)t , S respectively. 

6.3 Propagation of Discontinuity 

Proof of 2 of Theorem [2] 

Proof of ( |18[) : The proof is exactly same as in-flow case in Section 4.3. 

Proof of pob Recall that we have [h{to)]xo,vo 7^ foi' (^^oi^o) S 7o ^^'^ *o S iO,m.in{ti,{xQ, —VQ),ti,{xo,vo)}). The 

proof is exactly same as the proof of in-flow case in Section 4.3 except Step 2. We need to show a continuity of a 

boundary datum on 7_ U 79 . In bounce-back reflection boundary condition case, we need to show 

0= [ ft|[o,oo)x7_]to,:ro,i;o =lim sup \h{t' , y' , v') - h{t" , y" ,v")\. 

^^O t',t" e B(t;5) 

(a', ■"'). (a", ■"") e 7- n b({xo,vq); s)\{xq,vo) 

Because {y',v') is in the incoming boundary 7^, using the bounce-back boundary condition, we have h{t' ,y' ,v') = 
h{t',y' , —v'). Further due to the condition < io < th{xo, —vq) we have < t' < th{y' , —v') and 



h{t', y', v') = h{t', y', -v') = ho{y' + t'v', i;')e""^"'^*'"^°' -(Va^^)(-.1''+(*'--)"''-')'^- 

+ / {KJi + wT+i-, -)}(s, y' + it' - s)v', z;')e-''(^')(*'-*)-i'o' KyMS)(ry+(t'-r).',.')d-rf5^ 







w w 



and similar representation for h{t',y',v'). Using the continuity of v{^JjL^),Kyjh and wT^{^, ^) we have 

O=[h\[o^oo)x'yjto,xo,vo = Ijm sup \hoiy' + t' v' ,v') - hoiy" + t"v" ,v")\ 

*^o t',t"eB{t;5) 

{y'^""'), {y'\v") e 7_ n B{(xo, vq)-. S)\(xo, vo) 



xe 



-v(vo)to-Jg° v{^-^)(T,Xo + (to-T)vo,Vo)dT 



where we used the continuity of the initial datum Hq in the last equality. 
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